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ABSTRACT 


Real terrain is complex and its effects on the atmosphere are 
not well understood. The two-dimensional change-of-roughness problem 
is; the first step in learning, about the interaction, A critical review 
of relevant literature is presented, various analytical solutions are 
surveyed, and several numerical models are described. An experimental 
study of airflow from mustard, roughness length Zo=12 cm, to fallow, 
roughness z4,=1 cm, leads to several conclusions. The zero-plane 
displacement is found to be an important, yet neglected, parameter. The 
surface friction velocity overshoots and returns slowly to the equilib- 
rium value. The wind profile in the modified region is more nearly 
logarithmic than the Peterson model predicts. The Glushko model gives 
the best overall description of the flow transformation, although the 
Elliot theory, with the great advantage of simplicity, is in fair 
agreement with the data. Values for the ratio of the standard deviation 
of the vertical wind to the friction velocity (%,/us) reveal no variation 
with stability and are in good agreement with the value for laboratory 
flows. Values for the ratio of the standard deviation of the longitudinal 
wind component to the friction velocity (6y4/us) show a variation with 
wind direction, due probably to large-scale properties of the terrain. 
A plausible value is obtained for the constant of proportionality between 
the kinematic shear stress and the mean turbulent kinetic energy. A 
brief discussion of attempts to extend the theory to more realistic 


conditions is given, and some recommendations for future work are made. 
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CHAPTER I 
INTRODU CT ION 


1.1 The Problem of Inhomogeneous Terrain 

Micrometeorology has been uuch devoted to elucidating the 
fundamental relations which are presumed to exist between vertical 
fluxes of sensible heat, water vapour, and momentum in the surface 
boundary layer of the atmosphere, and the gradients cf the corresponding 
variables. To simplify the governing equations, the theory has assumed 
horizontal homogeneity, thereby employing essentially a one-dimensional 
model of the atmosphere near the ground. Implicit in the idealization 
is the existence of equilibrium between the vertical profiles and the 
underlying uniform surface, 

Natural terrain, however, is far from uniform. Apart from 
flat grassland, desert plains, or frozen lakes, few areas qualify as 
approximations to 'infinite planes'. Experimental sites suitable for 
testing the theory can be located only after extensive search. Most 
of the earth's surface is characterized by spatial changes in roughness, 
temperature, and moisture, as plainly evidenced by the chequered pattern 
Of agricultural lands. The effects of this 'patchiness’ on the vertical 
profiles and fluxes is yet to be determined, despite the obvious impor- 

tance of ascertaining the extent of deviations from the simple theoretical 
model. 

Because the fluid particles in the turbulent boundary layer 
possess a higher velocity parallel to the earth's surface than normal 
to it, the boundary layer is often said to have a 'memory'. Once 
equilibrium is established with a particular surface, the flow 
characteristics will tend to persist even after a new surface is 
encountered. Intuitively, it is expected that adjustment to the new 
surface will grow very much like a boundary layer on a flat plate. 

If point observations, taken from a tower for instance, are to be 
representative of these new surroundings, then clearly they must be 
made within the new equilibrium layer. Thus, it is of considerable 


practical concern to learn something of the response of the flow 
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properties] tostypicalssunface: irregularities..<In the! fieldia rule-of 
thumb in the form of a height to fetch ratio would be most valuable. 
Observation sites could then be placed far enough downwind from the 
beginning of any particular type of land cover to ensure that measure- 
ments up to any chosen height would be representative of that cover. 
For areas of very limited extent the maximum height of instrumentation 
could be calculated provided that the distance downwind from the start 
of the cover was known. 

Representation of the effect of the earth's surface on the 
atmosphere as a whole has long been an obstacle in numerical forecasting 
and general circulation modelling. For short-range forecasts it appears 
not to be important, but for medium and long-range forecasts it is vital. 
The surface is an important link in the energy exchange process; 
induced vertical velocities and the convective transfer of heat and 
moisture couple the boundary layer with the atmosphere above. While 
it is not expected that small crop to crop variations will be significant, 
the overall difference between land and sea, or mountain and plain, 
should exert some influence on large scale properties. Air mass mod- 
ification is one familiar manifestation. The formation of ideas about 
the effects of surface inhomogeneity is thus of some consequence on 


the synoptic scale as well. 


1.2 The Black Box Analogy 


While much of the interest in the problem of inhomogeneity 
has been practical, another motive of the investigation has been to 
gain a better understanding of the structure of turbulent boundary 
layers. Clauser (1956) has suggested that the turbulent boundary layer 
is a complex non-linear system resembling Maxwell's problem of a 
complicated machine in a closed room, connected to the exterior only 
by means of a number of ropes. By pulling the ropes and correlating 
the responses in an intelligent fashion, can the complete behaviour of 
the machine be predicted, say when a stranger pulls the ropes in any 
manner whatsoever? Like the machine in the closed room the actual 
mechanism of turbulence is inaccessible, the only connection to the 


‘outside' being a number of measureable parameters. 
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In modern terminology, Maxwell's room is called the 'black-box' 
problem. The black box has a number of inputs such that varying functions 
can arbitrarily be fed to the unknown mechanism inside. A series of 
outputs from the box indicate the behaviour of the mechanism inside. 

The objective is to establish what the output function will be for any 
given input. 

The turbulent boundary layer has as inputs various external 
influences such as pressure gradient, roughness, and heat flux, to 
which it is subjected in travelling along a surface. Typical response 
variables are the mean velocity profile, shear stress, heat transer, 
and turbulent energy. Recent investigation into the problem of 
inhomogeneous terrain is thus a search for a better understanding of 
the nature of turbulent flow by submitting it to a sudden change in 


one of its inputs and examining closely its response. 
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CHAPTER II 
THE INTERNAL BOUNDARY LAYER 


2.1 Formulation of the Problem 

The simplest situation involving a change of roughness is 
the two-dimensional case (Fig. 2-1), At x=0 there is an abrupt 
transition from a surface of uniform roughness z to a surface of a 


ol 


different uniform roughness z The surface whose roughness length 


o2F 
is the smaller will be called the 'smooth' surface; the one whose 
roughness length is larger, the 'rough' stirface, The first surface is 
assumed to be of sufficient extent upwind that equilibrium has been 
established between it and the various profiles. The stress will be 
effectively constant with height in the surface layer under consideration 
and equal to the surface value, Under neutral conditions the mean 
vertical velocity is zero and the mean horizontal velocity is given by 


the logarithmic law, The three characteristics of the upstream flow are: 


(1) cA Fae. pethateleemcnmarcons tant .=,_7, 


el 
D2 


where, t, is the kinematic shear stress, and) 2%, is the suriace value, 


(2) | u,= (u,, /k) ln 2/24 


1 is the surface roughness, Way is the 
oO ray 


Suriace frictiom velocity, k is von Karman’s constant, and z is the 


wiere wu. 1s the velocity, z 


vertical coordinate taken positive upwards, 


(3) w= 0 


where Wy is the vertical velocity. 

When the flow encounters the step change to a surface of 
different roughness, it is reasonable to assume that the flow responds 
to the new conditions, not at all levels, but only in a layer adjacent 
to the new surface. The height of this layer above the surface will 
increase downwind from the discontinuity in roughness just as a 
boundary layer grows on a flat plate, Because of this and the fact 


that a boundary layer already exists over the region, the new layer 


has been called the internal boundary layer. The problem is to determine 
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streamline through h 


Figure 2-1. Schematic diagram of internal boundary layer development 
for smooth to rough flow. 


the shear stress distribution, the horizontal velocity profile, -and the 
vertical velocity distribution, in the developing internal boundary 
layer, At large x distances the lower atmosphere (few tens of metres) 
should be completely in equilibrium with the new surface, and behave as 
if it had developed over an infinite uniform plane, 

The equations of mean motion for a two-dimensional, incom- 
pressible flow, in neutral stratification and steady-state conditions, 


neglecting the viscous and Coriolis terms, are: 
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Wherevp Ls"the pressure,” "is the density, u’ is the fluctuating or 
turbulent part of the horizontal wind, w' is the fluctuating or turbulent 
part of the vertical wind, In addition there is the equation of 
continuity, Pe ee ae, se 

2x oz 
which provides a system of three equations in six unknowns, To simplify 
the problem further the following assumptions are made: 

(1) The contribution to the momentum balance of the normal 
stress components ae and ws is small compared with other terms in the 
equations, Observations on the whole have confirmed this, © 

(2) The pressure is constant throughout the region of the 
roughness change. The atmosphere is, of course, always subject to at 
least a small synoptic pressure gradient, but over distances of 
micrometeorological interest, it will be negligible, Near the dis- 
continuity local pressure fields are expected, but provided attention 
is restricted to distances somewhat downstream the error should not be 
too great, This assumption will be re-examined later. 

(3) Initially w=0 everywhere, so the left side of the vertical 
momentum equation is zero, Because of uniformity and equilibrium upwind, 
the shear stress is not changing in the x direction so ke O. 

Hence the pressure term is also zero and the equation as a whole can be 
dropped from the system, Near the discontinuity, however, w is definitely 


non-zero and is changing rapidly in both the x and z dirictions. So too 


is the shear stress. Downstream from the step it is expected that the 
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usual boundary layer conditions will prevail, that is, all the terms in 


the equation will be small, At large distances the return to equilibrium 


will again set the whole equation to zero, Thus, except close to the 


discontinuity, the vertical momentum equation can be neglected. This 
assumption will also be re-examined later, 

Alternatively, the normal turbulent stress terms can be 
retained so that away from the roughness change the vertical equation 
of motion simplifies, after integration, to 


lah: 
P+ pw’ = constant = spe 


where the constant ie can be identified as the static pressure at the 
ground (Plate, 1971, p. 141), The préssure and normal stress terms in 


the horizontal equation of motion become 
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The assumption of isotropy, u' =w' , would reduce the right side to 


a single term expressing the fact that the pressure gradient is impressed 
by the conditions that determine ee In shear flow it is impossible to 
have isotropy (Stewart, 1956), but away from the discontinuity changes 

in pReiss with x will be small, Therefore the second term can be 
neglected, again revealing that the pressure gradient is controlled by 


conditions outside the boundary layer, namely, the synoptic gradient, 


With these assumptions the governing equations reduce to 
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a system of two equations in three unknowns, Solution of the set 


requires a third equation plus boundary conditions, 
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2,2 Qualitative Features of the Flow Modification 

When air passes from a smooth surface to a rough surface, a 
series of reactions are set up, A pressure field is produced and the 
air will be affected before it reaches the transition, The pressure 
will probably be felt only a few roughness lengths upstream, but, in 
addition to causing an initial deflection of the streamlines away from 
the roughness elements, it may cause local distortion of the velocity 
field, One example observed by Meroney (1968) is that near the edge of 
a €rop Or forest the velocities near the @round are greater than those 
higher up or farther inside the canopy. Plate (1971, p. 144) claimed 
that the pressure acts as a smoothing function to make the profile 
Changes less gbrupt, and that it is restricted to a fairly narrow resion 
near the discontinuity, He suggested that the overall effect may be to 
shift the origin of the internal boundary layer, much as occurs in jets 
and wakes where the boundary layer assumptions are not satisfied at the 
physical origin, and a 'virtual' origin appears. 

As the air rises tO “pass (over the roughness elements there 
may be a slight tendency for the fluid particles in the lower layers 
to accelerate, The increased roughness would quickly counteract this, 
the higher drag removing momentum, The velocity gradient would increase 
and cause an increase in the turbulent energy. In other words, some of 
the energy of the mean wind is converted into turbulence to satisfy the 
increase in momentum flux at the surface, The increased turbulence 
helps bring down, from successively deeper layers of air, the momentum 
required to overcome surface drag. The region near the discontinuity 
is one of intense turbulence production where a strong stress field is 
present. Meroney (1968) also observed that at the edge of a canopy the 
surface drag and turbulent intensities were increased considerably over 


the values farther inside, 
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2,3 Boundary Conditions 

The complete situation is shown in Fig, 2-1, The streamline 
displacement is represented by $ and the height to which the flow is 
affected by the new surface is denoted by h, Above z=h the flow is 
still characteristic of the upwind conditions, Near the ground, start- 
ing perhaps a little downwind from the discontinuity, a thin layer exists 
where the profiles are governed entirely by the local surface conditions, 
This equilibrium layer is denoted by h', Between h' and h there exists 
a blending region in which: 

(1) the velocity distribution changes gradually from the 
logarithmic one characteristic of the second roughness to that character- 
istic of the first roughness; 


2 2 
(2) the stress changes from Uys CO Uy,3 


(3) the vertical wind velocity is non-zero, 
Accounting for the streamline displacement, continuity of velocity at 
z=h requires that 


uy (h) = uy (h- §) 
= (u,,/k) ln (he §)/z, 


Below z=h the velocity distribution can be represented by a profile of 
the form 


u fu, = 1k) ln z/Z 0» + £, (z/h) 


2 


D is the friction velocity at the second surface, and f) is a 


function of the parameter z/h, the height divided by the new length 


where u,, 


scale, This 'blending' function must be such that: 
(LY Motus hh f= 0 
Below h' the profile is simply the logarithmic one characteristic of 


the second roughness; 
Cua tee.-an. £ (1) must fulfill the continuity condition 
uy (h) = (uy, /k) ln z/z + uy, £, (1) =u) (h-$ ) 
Similarly, the stress and vertical wind distributions can be written as: 


2 2 Z 


= - h 
T, ara Ux9 2 (U4 Uy) £,, (z/h) 
wp se £.,(z/h) 


where the functions f,, and . must satisfy the conditions: 
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(a) 2for gen. f,= 0 a 


The shear stress is constant with height and equal to the square of the 


friction velocity at the second surface, 


f Se = = 
ong We f, 0 Wo 0 
There is no vertical wind in the equilibrium layer, 
2 
b = = _ 
(ba paige £, (1) 1 tau, 
The shear stress equals the upstream value, 
= iy = = 
At 2=h f,( y= 0 w, 0 


There is no vertical wind in the upstream flow. 


Kinatly.. ate the eround, 2.2 


o2° 
2 
TT Tye Uyo 
uy = 0 
w= 0 


A complete solution of the internal boundary layer problem will provide 
the “distributions of u, (2) w (2); and t,(z) as functions of x, This 
will determine the growth of h and h' with x and specify the functions 


f\> fi and fo. 
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CHAPTER III 


ANALYTICAL SOLUTIONS 


3.1 A Semi-Empirical Approach 

The earliest investigation of flow over a change in roughness 
was carried out by Jacobs (1939) in a rectangular channel. Measure- 
ments of the mean velocity profile were made at various distances down- 
stream from a roughness discontinuity both for the smooth-to-rough and 
for the rough-to-smooth transitions. The shear stress distribution was 
then computed from these profiles by graphical integration of the 


combined equations of motion and continuity in the form: 


ee es) w ou = 25 = z 
ae + a sz = er ) at Fired!) 


(Some references to Jacob's work erroneously assume that he used the 
Prandtl mixing length technique, for example, Schlichting , 1968, p. 
615. In actual fact, Jacobs rejected this approach when preliminary 
mixing length calculations were found to be in serious disagreement with 
observations.) Jacobs concluded that the new equilibrium value of the 
wall shear stress is attained almost immediately, and that the change 
spreads upward with distance according to the exponential transition 
function eg(z/x) =a exp(-bz/x), where a and b are empirical constants. 
According to his observations,: a =1, b=11.6, and the resulting inter- 
polation formula for the shear stress distribution across a channel is 
e = Hee fr, + (%-%) exp (-11.6 2/x)] 
H 

A similar approach was taken by Lettau and Zebransky (1968) 
to interpolate wind profiles between an initial state over land and a 
final state over water. They proposed a Gaussian transition function 
of the form exp -a eine , whence 


2 
usu, + (u,~ u,) exp -a (z/h) 


Brooks (1961) interpreted Jacobs' results in terms of per 
cent transition and examined some of the implications. The internal 
boundary layer does not begin to grow immediately at the transition; 


the distance required for 95% establishment of the new surface equil- 
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ibrium shear stress is almost 100 barrier heights for rough -to-smooth 

flow and over 200 barrier heights for smooth-to-rough flow, Hence tree 

interference is to be expected for a hundred or so tree heights downwind, 

Brooks also compared Jacobs's formula with one by Jensen for the depth 

of a turbulent boundary layer developing from the beginning of a clearing 
a= ae rte 341) a : 

He concluded that for changes in shear stress of a factor of two a 

height-to-fetch ratio of 1/50 is sufficient for equilibriun, 

There is some question about the distances from the wall 
measured by Jacobs (Brooks, 1961) because the smooth surface was 0.3 cm 
heigher than the top of the preceding roughness elements to eliminate 
a pressure jump at the discontinuity, Also some doubt has been cast on 
the two-dimensionality of the flow, since the roughness change was 
applied only to the floor of the channel (Tani, 1968). Similar exper- 
iments in recent times by Logan and Jones (1963) for pipe flow, by Tani 
and Makita (1968) for channel flow, and by Antonia and Luxton (1971, 
1972) in a wind tunnel, all indicate that the wall shear stress does not 
attain its equilibrium value immediately, but that it overshoots and 
subsequently returns slowly to the equilibrium value, The velocity 
profile adjusts rapidly, establishing the law of the wall, The log- 
arithmic form, however, does not mean that the shear stress or the flow 
as a whole is in equilibrium with the underlying surface, 

Unfortunately, laboratory studies such as these are not directly 
applicable to the atmosphere, There are at least two difficult scaling 
problems: 

(1) The atmosphere possesses an overall, developed, planetary 
boundary layer, much thicker than the internal one formed by the roughness 
perturbation, In a wind tunnel, both boundary layers will continue to 
grow in the streamwise direction, The internal layer quickly grows to 
a considerable fraction of the outer layer; soon thereafter the tops of 
the two layers actually intersect, 

(2) The atmosphere has a thick frictionless region and no 
rigid upper boundary, The wind tunnel has a fixed upper wall and a 
rapidly diminishing 'free flow' region, Also the fixed volumetric flow 
rate allows only the shape of the profile to change whereas in the 


atmosphere the whole flow can accelerate (Brooks, 1961). 
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3.2 A Momentum Integral Technique 

The horizontal momentum equation together with the continuity 
equation can be integrated across the depth of the internal boundary 
layer after the manner of von Karman to yield the integral momentum 
equation hk 


h 
ebonfdee Spd ~% 2 f ude Ziut Mast ste 4, A = 
dx dx 

az 


2 


This expression simply states that the horizontal flux out of a region 
bounded by the top of the boundary layer, the surface and two vertical 
lines a distance dx apart is equal to the net gain of momentum due to 
vertical flux into the region. The latter is given by the difference 
in shearing stress between the top and bottom of the column considered. 
For closure Elliot (1958), in a pioneering paper on the 
subject, assumed that the velocity profile for z<h immediately takes 


on the logarithmic form of the new surface roughness, 


u, = (u,,/k) ln z/z 0, 


implying that £ (a/b) = O. The boundary conditions were those previously 
described except that velocity continuity was imposed in the form 

u,(h) =u, (h) at z=h, tacitly assuming that the streamline displacement 
is negligible compared to h. The formal solution is extremely clumsy, 
but Elliot found that the growth of h with distance is well represented 


by 


Sats) ) 


forex / 2a >3L0 


) 02 


h/zissa Ge/z 


o2 o2 


Near the origin the method could not be expected to yield accurate 
results owing to pressure effects, The coefficient a is a slowly 
varying function of the roughness change parameter M=In(z 04/259) 
which can be approximated by 
a= OD. Joed D.03.i 

Hence the height of the internal boundary layer is greater for air 
flowing from a rough-to-smooth surface than for the reverse, The slope 
of the interface is about 1/10 for distances of micrometeorological 
interest. 

A surprising feature of this growth equation is that it does 


not contain either the wind speed or the friction velocity. In 
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boundary layers growing on a flat plate the height is dependent on some 
power of the free stream velocity. For turbulent flow over an aero- 
dynamically smooth plate the growth equation is given by 


0.8 -0,2 


hax (ase ) (Schlichting, 1968, p, 399). 


and the velocity profile is u/u, «In (u,z/0), For flow over a rough 
plate the velocity profile takes the form u/u,%ln (z/z,). A 


comparison reveals that the surface characteristics measured by z 


w 


have replaced the viscosity factor ius. as the relevant parameter, 
If this can be extended to the growth equation and wv is replaced by u,z 
then 
OLS -0,2 
h Ce xX (u, /4,2,) 0, 


0 0,8 


and (u, /u,) z h/z OL (x/z 0) 


B 
ut u, fu, « ln (h/z.) 


On 0.8 


so (ln h/z) h/z. x (x/z,) 


2 
The factor (ln ae is virtually constant over a large range of 
h/z Varyvine irom 1.4/7 At B/Z = 10" to 1.79" at icy seit. Therefore 


h/z a (x/2,) 2 


This is the dimensional argument put forward by Elliot in support of 
his growth equation, 

The Elliot theory seems to incorporate the essential features 
of the flow problem as it gives results in general agreement with 
observations (Echols, 1972; Panofsky and Petersen, 1972). Also see 
Fig. 3-2. His basic assumption, the immediate attainment of the 
logarithmic profile, has been demonstrated experimentally both in the 
laboratory, as already mentioned, and in the atmosphere, by Bradley 
(1968). Nevertheless, the solution is not complete because the shear 
stress distribution in the vertical has not been specified. To do this 
a relationship between the shear stress and the velocity profile must 
be assumed. If Jacobs' work is any indication, a mixing length assump- 
tion will not be adequate. The elegance of Elliot's simple approach 


becomes most evident as other theories are examined. 
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In an attempt to improve upon Elliot's model, Plate and Hidy 
(1967) assumed that the streamline displacement is not zero, but small 
compared to h. The continuity condition at the interface was then 


written as 


UW /Y yy = (i/k) (ln h/* 4 = te fat) 


by expanding the logarithm in a series and dropping second and higher 
order terms. A simplification was obtained by integrating the momentum 
equation between x= 0 and x=x, the shear stress being set equal to an 
average value along x, and the streamline through h(x) assumed to be a 
straight line. Their results are incorporated into a nomogram (Fig, 
3-1) from which h/zo4 and we 


x2 
limited to short fetches where the blending region is not a significant 


Can be obtainede | Lars modvrrteat1tom as 


portion of the internal boundary layer. 


3.3 A Simple Diffusion Theory 
The experimental findings of Dyer and Pruitt (1962) prompted 


Dyer (1962) to assess the rate of adjustment of profiles and fluxes 
using simple diffusion theory. For near neutral stability Philip (1959) 


found a solution of the diffusion equation 


u se a¢ $ = —- de 
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where © is the concentration of the diffusing entity, ® is the 
corresponding vertical flux, and K is the eddy diffusivity. 


Assuming that the wind profile is given by the power law 


usu (2/2)! 


where u. is the velocity at a reference height z taken to be 1 cn, 
r 


‘and assuming that K is a function of z and not x given by 


6 
K=kK. (2/25) /7 
r ic 
the equation became 
pag 2 2(t 35) 
Te Ox = 22 Y Bz 


Solutions were found in terms of the parameter 
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(a) Internal boundary- 
layer thickness 


(b) Ground shear-stress 
distributions 


Nomogram for calculating internal boundary layers (from 


Figure 3-1. 
Plate and Hidy, 1967). 
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F = (49/81) (u MK) (2/x) 7/7 


The rate of adjustment was determined by considering the ratio / Qo 


where .. is the flux at x=, that is, when full equilibrium has been 


established. Thus, 


b/ta = 25 / (2). 


The profiles experience an adjustment rate equal to that of the fluxes, 
which follows from assuming that K is independent of x. Philip's 


solution gave 


fe tee ei 


where I(f ,-1/9) is a form of the incomplete gamma function, the parameter 
-1/9 arising from the power law assumptions. The solution is a function 
of x, 2, and u/K ,» but Dyer noted that the ratio u /K is fairly 
constant over a wide range of wind speeds. Using Philip's tabulation 

for the solution the height to fetch ratio for 90% adjustment (0/02=<9) 
can be computed (Table 3-1). 

Dyer's intention was to provide working estimates for research 
into heat and moisture transfer, as in the later work of Rider, Philip 
and Bradley (1963), and Dyer and Crawford (1965). Although specifically 
excluding a change of roughness in the original treatment, Dyer (1965) 
considered it to be an approximation to shear stress and wind profile 
modification. It is difficult to accept this view because his theory 
assumes that K is a function of height only. By definition K= ku,z 
and in the change of roughness problem u,.varies with x; it follows that 
K mst also vary with x. The adequacy of such a simplified diffusion 


equation is also questionable. 


: 
ho. Sid | 

= ey Raa aeaten 4 

a alpen ot aiteablenas ed toniiirretl sa teens hf ty 10 & 

seed 064 oi Oh Lian fh) aay eh, at te Is xulh: ether 


, tote - «Nad 


; 7 = a 


aciy\ s&s 
@ (23) \ a ~, 5 > ; 7 


Cm s 


a 
peeeut? ate le Jali ou? fp4ne.e0n? s4unteul la ph quo lrape aot rong ie 
7 


e'yliant ot Te GHesHISSAL GL VV ONT Pieesh Geter em aie dab a 
: EE wat ? te 


j T H A a : 
® eu - — 7 
on hinte oq ald metsanir anvitg Au imcoys Gila FO te? ' a St | pM wth 
miainwt i at rv’ 1isY ab of pauies@ mol eanty Gers Mayer priaaee @\ cme 7 
yfoiet @F Ay: Udand i) Jvita fay! Pay Le, oats tine a ye « 
mesniitad @'2i0909 girl higwirh tnky 1 oe; iy -& tHYu sank 
(*. me TT ?La SO ae is iis j¢ 9 {hi i) viel dettin wa a 
(iit) pete « vie 7 
Avaswenss 100) Ca eeeesedg pierce Ghivecd us Gale ota emia o! 
i pat “nik Procw [ sale 4 ~ ,102 y ieton | swudl 64 
vtlestiicoce Gyyotiin .COAMT) OSretwiyo/Tiw vey? Dan wer) orl = , 
(202!) 36 taguwiethr? lapighea aii ni * ery 1S, GG + inthe 
af 456%) Dniw Gib Sas72a° 1140 270164 Bi ere A je So ey We week hewn: 
eas eft C8ay py era? Wss16 45 (fumes VXEe « c 371i .,. vp enelIes 
t,08 2 noe dbada ib  pwteo thglot jy ake reel be 4 gents e ane 7 
sab evelind 4% iv ftir SAl1Se@V <i cides sitelion Ta aunt af 2 


aes Matti lquke 8 dean to YOnvgen. oft” Hiyfdet vise onde amet © 
. ' ,didenetseuep cate a wake ix : 


18 


Table 3-1, Values of fetch and height for 90% adjustment. 


Height (m) Fetch (m) Height/Fetch 
Oxo 70 1/140 
L.0 170 1/170 
pid) 420 1/120 
Sell Lb; 350 1/270 
LOO 3, 900 1/330 
20.0 8,100 1/405 
500 Z2Ona00 17/530 


ol col 


Figure 3-2. Schematic representation of the transition region. 
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3.4 An Advection Theory 

In a dynamical approach R, J. Taylor (1962) identified the 
region of stress modification with that of streamline displacement as 
shown in Fig. 3-2. Here AC and BD, both lines of constant stress, mark 
the boundaries of the transition region, and CD is part of a streamline 
in the rough-to-smooth change. A streamfunction ¥ can be defined so 
that the horizontal momentum equation becomes 


a oe eet ees ORL) 


we rc) (a) Se. One Oz 


Integration over the region ABCD gives 
Ww 


f(g a) d% =f 4 dx 


ABCD 
where usu) is taken along the streamline. 
Accepting Jacob's results, Taylor then assumed that the distance between 
A and B is mich less than the distance between C and D, making ABCD 
approximately triangular. Hence, 
— eg es pe St meee. 
ABCD 


where X is some x-coordinate between C and D (about midway if the 


shear stress varies regularly along CD), and 
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ce) 


For neutral conditions the stream function will be 
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u and z , the parameter “X can be evaluated 
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by substitution for the integral. 


Thus, given u 


In practice, however, only three of the four are known and 


Uxo is to be determined. Taylor assumed that 
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and together with the previous equation obtained an expression for x 


in terms of the roughness-change ratio m=exp M = Zo /% 0° 
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If AC is taken to be vertical and the streamline displacement 
is assumed symmetrical about x =X , the new velocity profile should be 
fully established to a height z at a downwind distance of twice X(z). 
For typical values of m this leads to height/fetch ratios of 1/100 to 
ITSO. 

Taylor's own wind tunnel experiments revealed some variation 
in the surface shear stress after the transition, but equilibrium was 
quickly reached, in confirmation of Jacobs' results. Aliso the stream- 
lines were observed to drop off at or a little before the transition, 
indicating that AC was indeed perpendicular. However, as Panofsky and 
Townsend (1964) pointed out, almost all the measurements were made in 
a region where the internal boundary layer actually intersected the 
tunnel boundary layer. 

This solution is not, of course, a complete one; no information 
is supplied about the profiles in the blending region. The parameter X 
is the scale of the transition and serves only to provide practical 


estimates for field work. 


3.5 A Parabolic Shear Stress Distribution 

Under the mixing length assumption, namely, 2S ihe! 2 ae i 
the logarithmic profile used by Elliot in the internal boundary layer 
implies constant stress es up to z=h and then a discontinuous jump 
to ie To remove this unappealing proposition Panofsky and Townsend 


(1964) assumed that 2% is a linear function of height such that 


oe == x9 2 (uyy os uy.) (Z/d) 


where d is the depth of the internal boundary layer. 


With the mixing length relation 
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the downstream profile is thus 


u,/u,. = C/k) ln Z/Z oo + (U4 yg) /k (z/d) 


the second term on the right being the blending function £(z/d). 
The conditions of incompressible, steady flow, and velocity continuity 


at z=d, together with the conservation of momentum, yield an ordinary 
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differential equation for the growth of d with x. For values of d that 


are not too small the solution is given by 
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where x 6/2, 20 foriM= 3etand forms 


p 
aka Jz, = exp (14M) (M/18 -64/27+88/27M + 67/9M7) 


One consequence of this continuous stress assumption is that 
the definition of the interface between the two flow regions is changed. 
In Elliot's theory h is the intersection of the two logarithmic profiles, 
whereas here d is the top of a transition zone, and so is somewhat 
higher. Except for short fetches (where neither theory is valid anyway) 
the theories are in close agreement. Near the top of the internal 
boundary layer the profile in the Panofsky-Townsend theory begins to 
deviate from the lower logarithmic profile and enters the upper logar- 
ithmic profile smoothly. The height/fetch ratio for 90% adjustment of 
shear stress is thus somewhat less than in the Elliot theory and is 
about 1/20 (Townsend, 1965; Dyer, 1965). 

In assessing the assumption of the equilibrium relation 
v= kz(du/dz) Panofsky and Townsend pointed out that the time necessary 
for substantial adjustment of the Reynold's stress is of the order of 
the turbulent energy divided by the rate of energy production. The 
distance travelled while the adjustment is made is nearly 


x = Ea /(u, du/dz) where E Wbte + aye + le 


The surface separating adjusted from unadjusted fluid is 
x/z_=(E/u,”)(d/z_) In d/z 
fo) * oO fe) 


which has much the same form as the previous solution. 
Another line of argument is that of Miyake (Panofsky and 
Townsend, 1964) who assumed that the vertical rate of spread of the 


surface change is proportional to the standard deviation of the vertical 
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wind, d(d/z 9) /d(x/z ) =A6,/u . With the assumptions 6, = uy, 
and u=(u,/k)1n d/z this yields 


x/z = CL/Ak)(d/z 0) (In d/z - 1) 


Again the solution has much the same form as that outlined by Panofsky 
and Townsend. 

Observations by Blackadar, Panofsky, Glass and Boogaard (1967) 
of the interface height agreed reasonably well with both the Elliot and 
Panofsky-Townsend theories. The extended validity of the relation 
~"*=kz du/dz in near neutral but nonequilibrium conditions was also 
investigated. The non-dimensional wind shear b=(kz/ 2+") du/az was 
found to deviate from the value 1.0 in a nearly linear fashion with 
height, the error at 1/10the interface height being 5%. Measurements 
of the surface shear stress on a beach taken by Hsu (1971) were also in 
fair agreement with the predictions of the Panofsky-Townsend theory. 

A series of atmospheric experiments related to the change of 
roughness problem has been performed with artificial roughness elements 
over the ice of Lake Mendota, Wisconsin. Kutzbach (1961) used a 
triangular array of bushel baskets, which may not have simulated a truly 
two-dimensional situation. His results were in good agreement with both 
the Elliot and Panofsky-Townsend theories (Fig. 3-3). Stearns and 
Lettau (1963) and Stearns (1964) used Christmas trees in a circular array, 
two-dimensionality again being questionable. These experiments were not 
designed for the purpose of testing the internal boundary layer theories, 


so it is not surprising that the data are inconclusive. 
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Figure 3-3. 


Comparison of bushel basket data with the theories of 
Elliot and Panofsky-Townsend (from Panofsky and Townsend, 
1964). 
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3.6 A Similarity Theory 

To improve the estimate of the blending function fy Townsend 
(1965, 1966) investigated the conditions under which it can be assumed 
to be a similarity profile dependent only on a length scale 1: He 
assumed that the velocity downstream of the discontinuity consists of 
three parts: 


(1) the profile u. upstream of the discontinuity; 


(2) a component he that has to be added because of streamline displace- 
ment ; 
(3) ‘the remainine ditference u’ due to acceleration of the fluid in the 
internal boundary layer. 
For streamline displacement &(z) small compared to 2Z, 
Similarity entered 


2) Sal, % /kz from a Taylor series expansion of u 


L 
through the assumptions 


' 
a= (uric). ECs fe) 
s s 


2 Ys Z 
aM RC ORC IS Shy) KC 2/14) 


where u. is the velocity scale. The scales 1. and ie depend only on x. 
The functions f and F were assumed to depend on z/1. and NOE Jon x. 
They were related through the mixing length so that, putting Nees 
cL) ery, 
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The equations of motion and continuity were combined to yield 
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The two scales were then found from the conditions that f and F are 
independent of x (in the sense that they depend only on 1), and that 
the velocity profile near the surface is logarithmic. Townsend's 
original work contained a small error in not satisfying the latter 
condition. Blom and Wartena (1969) corrected the derivation to remove 


this inconsistency. The revised scales are given by 
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+i 
a 
U = “Uy M/(14P) [ in Usions, Teeter (ee bd 


-1 
where P=(1n 1i/z 4 -1-y) and Y= Euler's constant = 0/577 


The downstream surface shear stress can be calculated from 


and the velocity profile is given by 


u,=Uy + Cu /k) [Bi(-1 )C1+P) - P(1 - exp oe 


While 1. is the same order of magnitude as h, the adapted layer according 
to Blom and Wartena is of order 0.1 1s This results in rather long 
adjustment distances, the height/fetch ratio being about 1/300. Fully 
developed boundary layers, even to heights of 10 m, will be difficult 


to find anywhere if this assessment is correct. 


3.7 An Exchange Coefficient Theory 
Russian work on the problem of inhomogeneous terrain is. 
based on the exchange coefficient K=ku,z (Panchev, Donev, and Godev, 


1971). The horizontal momentum equation in the form 
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If the perturbation in wind velocity is small compared to the equilibrium 
velocity over the second surface, linearization is possible, and an 
approximate solution can be found by the method of Shwetz (described in 
Panchev, Donev, and Godev, 1971). Gandin did this for the case of a 
sudden change in surface friction velocity, but not roughness. 
Donev was able to solve the non-linear equation of the general 
case also using the method of Shwetz. The solution for the blending 
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and the internal boundary layer height can be determined from 
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The solution of the linearized equation differs only in the absence of 


the last terms in both equations. 


3.8 Experimental Evidence 


The best atmospheric data presently available are those of 
Bradley (1968) who measured simultaneously both velocity profiles and 
surface shear stresses. His experiments involved flow over tarmac 
(z= 0.002 em) and an artificial roughness mat made of vertical spikes 
(20> 0.25 cm). The surface shear stress was measured with a drag plate. 
The results are shown in Figs. 3-4, 3-5, and 3-6. 

Several important conclusions can be reached. The velocity 
profile is practically unchanged above the point where the lower profile 
intersects the upstream profile, implying that the streamline displace- 
ment is negligible. Also there is little deviation from the logarithmic 
profile over most of the internal boundary layer. Near the top of the 
layer the profiles do exhibit a slight curvature. The discontinuity 
is rather sharper and the top of the modified region is iower than that 
predicted by the Panofsky-Townsend theory. The 4/5 power law is 
confirmed and the predicted assymetry between smooth-to-rough and rough- 
to-smooth flows is not apparent. When the top of the internal boundary 
layer is taken as the upper limit of the linear portion of the profile 
the slope of the interface is about 1/20. In fluid mechanics the 
‘law of the wall' is usually assumed to extend over the lowest 10% of 
the boundary layer. With this in mind Bradley recommended that the 
height/fetch ratio for the equilibrium layer is 1/200. 

Most of the stress change occurs within a short distance of 
the transition, more rapidly than either the Elliot or Panofsky-Townsend 
theory predicts. The recovery from the initial overshoot is less rapid 
in the rough-to-smooth transition. There is a hint of a minimum in the 


smooth-to-rough case. 


= 


a 7 
x 


; ar a oe sida weynt wan beinnd | ‘6 
2g «fee (i. ahaa eo 


i id geo 


Cyae yo 1 : 7 j A e _ 7 
Z' 7 ° a i] a = 
\ a oh 4, 6 tis fa | : S) 
= : : 
a J ee Vile ’ ;. &s4 = [ide =) a> 
i.e rae - oo vy ie \ o® 4 -_ 
» cooneeds 249 i YEA vis pabeeap]e Het ean One le wlites 
Are Gilips AO Nd aT ioat Poe 
: 7 
_ 
ove is? (nteitzes SE: 
mala ates 0 aid. sak aA? : 
i aye i? ears Ff ; ove 1 50G8) eolber 
eee 
: r : cy am Le bu j (je v9 = 1.0 = - 
£ Pa oe a ' jh. cae been it «(ae Oz s 
lieth Bilin 6 eet nl eas peels eee 
kes wi NSIT U7 PS OT eT carte ae im: leaveuee > . 
ph re Sail a) Tua | austere oi PQe Coote et ef(lesq 7 
‘ ; > —— 7 —_ 
é ( ie tm » a Gath @ ptramy 9M 20 yen toy eh) €29eeTe 
= 7 . 
— pe Poe | if) ae) oes sgt is;) i Le] fem ivi 2 aa as 7 
ah 16 get sil reg 4A aad lerreosgl orl Se. 1 om Vee otsio aq 
eect Most bh att ieevews Piside w id ot rh se ae royal Z 
+5 . 
judy sel t Thar | ‘J rm 1 ot j ig] GY Sta. © yale vatteat ai 
al wal pagtyey ia sei VOTH SAVE fees wiowdé't tha ye eedalt 


i 
‘ a 
= 


, 
~ ener? bine iment-©) 9008) imGidaalh yer with na aeEr es @it® hese oman 18 
yeeenesl lanvetai wilt fo 7O! qs ag, anegqa hor wd wie os Genie ste 


' 
La bk 4 Ye cmelirwon ss VI. aly bis Teint! 1 O40 y. oun ~ oNton - meal 
if@ mys quiere ask ni eve a » santos tat Te i , lis 


a 
to TD) saqeed ait) teyy Imngrrs yo) homemaa YL Ines. a) ite ity wd 


wkd vats. oWiOmOaet gaThiia Vata «at: sible Agiy send: assed ut 
ee oud tl Sre if? a0) oun vada 


to asp sa thcaveds a pad lw eiui & nisi ceria yt 


7g] 


Wie h (6.42) P-H 
2 oe! 


Wages 0.038 tos, 


0.4 0.6 0.8 1.0 1.2 
U/U529 
Fetch over spikes, m 
+ 16.42 
4 6.42 A 16.42 
Data 4° 2.32 P-T theory, B 6.42 
2 nik} Ca 2332 
e 0.32 ; 


P—T = Panofsky—Townsend (1964) 
O P-H = Plate—Hidy (1967) 


u2>/u2, 


P = Peterson (1969) 


0 4 8 12 16 
Fetch over rough surface, m 


Figure 3-4. Experimental data of Bradley (1968) compared with various 
theories for the smooth-to-~rough transition (from Plate, 
1971, p. 142). 
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Figure 3-5. Experimental data of Bradley (1968) compared with various 


theories for the rough-to-smooth transition (from Plate, 
LO7ie pe Lo). 
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CHAPTER IV 


NUMERICAL SOLUTIONS 


4.1 A Mixing Length Model 

Rather than assume a special form for the shear stress or 
velocity distribution, P. A. Taylor (1969a) chose to solve the govern- 
ing partial differential equations by numerical means, using the 
mixing length relation for closure. To facilitate the computations, 
slightly different co-ordinates were employed, the origin of the z-axis 
being shifted to the local roughness length Zoi" The mixing length 
then took the form Q=k(2tz..), the initial profile became 
u/uy., = C1/k) ln (2tz,)/251> and the boundary conditions were given by 
u=w=0 at z=0. Velocities were scaled with respect to u,, and lengths 


with respect to Z 9° Eliminating ~ the equations of motion were written 


as Ht Ni ee eee! 
OX avy ae Noes 
che +e re Se 
dx xy 
= = ss = =F = 
where U u/Uy4> W w/Uy4> bs x/Z 53 and. f= In. (z Z o9)/Zo9 


Setting 2= Z/Z 9 the initial conditions became 


U=Cl/k). In (22 4/204 + 1) 


Replacing the ¥ derivatives by finite differences produced a system of 
ordinary differential equations for dU/dX, which were solved by the 
Runge-Kutta method. 

The solutions were compared with the Elliot, Panofsky-Townsend, 


and Townsend theories as well as with a solution based on the form 
glee Uns + (ye) “Uyey) [10¢2/h)>- Is(c/ nee pty a) | 


derived by Taylor using the Pohlhausen technique (Goldstein, 1938, p. 
158) of choosing the highest order polynomial that will permit reduction 
to a single ordinary differential equation. The results of the numer- 
ical computation were in close agreement with the Panofsky-Townsend 
theory, especially the velocity profiles (Fig. 4-1). The shear stress 
profiles were also in good agreement with the self- preservation theory 


of Townsend (Fig. 4-2). 
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Figure 4-2. 


z/2o2 


20 


15000 


10000 


5000 


2 ee Panofsky-Townsend 
Sr ae self-preservation 
= 322)" polynonial 


en mixing length 


Comparison of shear stress profiles from mixing-length 
model with various other solutions (from Taylor, 1969a). 
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4,2 A Turbulent Energy Model 

One weakness of all the foregoing theories is the assumption 
of the equilibrium mixing length relation between shear stress and the 
velocity profile. The AFOSR Boundary Layer Conference (Kline et al, 
1968, p. 425) recommended that the shear stress be 'unhooked' from the 
mean velocity profile, in other words, that there is no unique relation 
between the velocity profile and the shear stress. Plate (1971, p. 155) 
stated that the mixing length is a purely local quantity whereas the 
shear stress develops over the whole of the boundary layer, and only 
its rate of change is determined by local conditions. 

To allow for this 'history' effect and to avoid specifying 
a stress-profile relation, Peterson (1969a) closed the system with the 
turbulent energy equation, using a number of assumptions to relate the 
various terms in that equation to the shear stress. The turbulent 


energy equation neglecting the pressure diffusion term is 
u dE pee one = Ou _ & oe. + 
Dx (ar4 ch sale ) € 
where E is the mean specific turbulent kinetic energy, E' is the 
turbulent energy and € is the dissipation rate. The following 
assumptions were made: 


(1) The shear stress is proportional to the kinetic energy, 


to eee ee 
where a in an empirical constant, The value was taken to be 0,16, 
(2) The dissipation is related to the shear stress as in a 
logarithmic equilibrium layer where production equals dissipation, 


aan) ee =vlé 


(3) The energy diffusion is of the gradient type, 
w'E! = -K,dE/dz 


(4) The energy diffusion coefficient is the same as that for 


momentum, 


K,= a t /(2u/ez) 
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The turbulent energy equation then became 


ry at / oe Ou 3/ 
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The boundary condition on the shear stress was imposed in the form 
2 
Ais [kz (du/dz) | 
© o Zs 


The numerical solution revealed several important features. 
The 4/5 power law for internal boundary layer growth was found to hold 
with the height h defined as the point where the stress reached 0.1% 
of the upstream value (this ‘was about equal to the height of the inter- 
face in the Panofsky-Townsend theory). Such a definition for h meant 
that the velocity was essentially the upstream value above z/h=0.8. 

The slope of the interface for large z/Z 4 was 1/10 and the Townsend 
prediction of a self-preserving shear stress distribution appeared 
valid to first order. The stress was closer to the upstream value in 
the upper 75% of the transition layer. The stress and turbulent energy 
adjusted rather more slowly than the velocity profile; only about the 
bottom 10% of the layer was in equilibrium. Hence Peterson suggested 

a neieit/ fetch ratio of 17100; 

There were two significant differences between the predictions 
of the Peterson model and the other theories. Firstly, the velocity 
profile contained an inflection point just below the interface (Fig. 4-5). 
There was a suggestion of this in Bradley's data (Fig. 3-4 and Fig. 3-5). 
Secondly, the non-dimensional shear 0 = kz/2!h( u/2z) deviated consider- 
ably from unity in the transition region. For the smooth to rough 
change >1, and for the rough to smooth change O<1 (Fig. 4-3). 
Observations by Busch and Panofsky (1968), and Yeh (Peterson, 1969b) 
support this conclusion. Deviations from unity were also noted by 
Blackadar et al (1967), as mentioned in Chapter ilI. Peterson asserted 
that this was because of the slow adjustment of the turbulent energy 
to the new boundary conditions relative to the wind shear, and the 
sharp change in slope occurring in the wind profile just below the 
interface. An important consequence is the fallacy of exchange coef- 


ficient or mixing-length theories. Because K=ku,z/ and ) Al as assumed 
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Figure 4-3. 
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Non-dimensional wind shear in the internal boundary 
layer (from Peterson, 1969a). 
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by these theories, calculations based on them could lead to serious 
errors as Peterson (1971) has demonstrated. 

Peterson's basic assumption is actually an old one, advanced 
by Townsend (1956) and recently by Bradshaw et al (1967). Bradshaw 
determined the constant of proportionality to be 0.30 (which Peterson 
mistook Fors 0/215) Lumley and Panofsky (1964, p. 128) obtained 
values of 0.12 from Brookhaven data and 0.22 from O'Neill data. 

Values supplied by Hinze (1959) are in good agreement with the 0.16 
used by Peterson. 

As for the second assumption, Monin (1959) has argued that 
the relation between dissipation and shear stress is valid even in 
non-equilibrium conditions. Pasquill (1972) has affirmed it on empirical 
grounds. Bradshaw (1967) wrote the relationship in the form € = mete 
where Le is the dissipation length, and noted that the results depend 
significantly on the form assumed for Le. Little is known about the 
form of L. for atmospheric turbulent flow (Shir, 1972). When production 
and ‘dissipation of ‘turbulent ‘energy are equal le = 2 =kz 
The wind ‘tunnel studies»of Antonia ‘and (Luxton (1971, 1972) have: indicated 
that both the dissipation length and the mixing length are altered from 
the value kz in the region near the roughness discontinuity (increased 
for rough-to-smooth, and reduced for smooth-to-rough). 

The last two assumptions, Peterson has argued, are not 
important since the major effect of the energy divergence term is to 
smooth the velocity profile without modifying its basic shape. Various 
forms were tried with no significant differences in results. P. A. 
Taylor (1972) has shown that the assumption Kak(zte et? simplifies 
computation. The results of this modified Peterson model are very 
close to those of the original, Details of a numerical-method are 
presented in Appendix B, 

Peterson (1972a, 1972b) has stressed the importance of dif- 
fusion of turbulent energy in the change-of-roughness problem, a feature 
of the flow described by the turbulent energy equation. Antonia and 
Luxton (1971, 1972) also found in their wind tunnel experiments that 
turbulent energy diffusion played an important role in the flow modifi- 
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4.3 A Vorticity Model 

The fact that Peterson's results did not agree very well with 
Bradley's observations for rough-to-smooth flow (Fig. 3-5) led Shir (1972) 
to include the pressure-gradient term and the vertical equation of 
motion. Neglecting the viscous and normal stress terms as before, Shir 


combined the two equations of motion to yield the mean vorticity equation 


where “ is the two-dimensional vorticity. The pressure gradient does 
not appear because the density variations are assumed to be small. 
The second equation, derived from continuity, was the streamfunction 


equation 5 2 
ihe ate a = Ww 
oar azi* 


where ¥ is the streamfunction. 
The third equation of the model was the turbulent energy equation as 
used by Peterson. Following Peterson it was assumed that: 

(1) ts ak 
where a was taken to be 0.22 on the recommendation of Cramer (Shir, 1972). 
Shir mentioned that Harsha and Lee obtained a value of 0.3 which was 
also used by Bradshaw (1967). Bradshaw and Ferris (1968) believed that 
a is a slowly varying function of the ratio of production to dissipation 
with a value of 0.3 in laboratory flows, but possibly smaller in the 
atmosphere. 


(2) a= poe Aa 
with the dissipation length k, given by 


ih = Piles atavy) ae 
& Dt 2 ee 
p+ 2 (2-2,)/> 


i > 2, 


where Za is the height of the surface boundary layer, taken as 10 m, 


and >» is the mixing length at the outer edge of the boundary layer, 


taken to be 40 m. This is a modified form of Blackadar's (1962) formula 
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for mixing length. Bradshaw found empirically that Le= kz near the walt 
and tends to a constant away from the wall. For atmospheric flow there 


are insufficient data to determine the exact variation with height. 
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The energy flux was assumed to be of the gradient type. Once again the 
results were found not to be sensitive to the precise form of the flux. 


The boundary conditions imposed were: 
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The results of the numerical computation were in many ways 
similar to Peterson's. The velocity profile had an inflection point, 
and the non-dimensional shear differed from unity in the transition zone. 
There were also several differences. For the rough-to-smooth change, 
the calculated stress values were a much better fit to Bradley's data 
(Fig. 4-4a ). For the smooth-to-rough case, a local minimum appeared 
in the downstream surface stress (Fig. 4-4b). The latter was intimated 
by Bradley (1968) in the discussion of his experimental WORK; but not 
pursued because no theory at that time predicted it. The velocity 
profiles, however, did not fit nearly so well (Figs. 4-4c, 4-4d), and 
led Shir to question some of Bradley's data. 

Shir also pointed to the inconsistency among authors in the 
definition of the internal boundary layer height h. Some chose to 
define it in terms of specified deviation from the upstream velocity; 
some, in terms of stress. Shir did both, defining a stress boundary 
layer height as the height where the stress is within 1% of the upstream 
value, and a velocity boundary layer height where the velocity is 
within 1% of the upstream value. The latter compared favorably with 
that used by Bradley, and the 4/5 power law was again upheld. The 
height of the interface was found to increase with roughness ratio m, 


but the slope remained the same. This suggested a relation of the form 
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smooth flow 


Figure 4-4, 


ShareGtrom Shar, 1972). 
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The slope was larger for the smooth-to-rough transition. The height of 
the stress layer was about twice the height of the velocity layer, 
implying that the velocity adjusted more slowly than the stress. The 
height/fetch ratio for the stress layer was 1/10 while for the velocity 
layer it was 1/20. 

The shear stress in the transition region was closer to the 
upstream value for the upper 75% of the stress layer, as pointed out 
by Peterson. However, if the velocity height was- used, this was 
found to be true only for the upper 50%. The shear stress distribution 
was not linear, nor was it similar at different downstream distances. 
In the smooth-to-rough transition the shear stress was constant with 
height (to within 10%) for the lowest 1/10 of the layer. For the rough- 
to-smooth transition this was true only for the lowest 1/20 of the layer. 
Because the velocity profiles were logarithmic for the lowest 10% as 
well, this was taken as an indication of equilibrium. The height/fetch 
ratio was thus 1/100 for the smooth-to-rough case and 1/200 for the 
rough-to-smooth case, away from the discontinuity. 

Peterson (1972a) has emphasized the similarities between his 
approach and Shir's, adding that it is difficult to judge the merits 
of a model using just one set of observations. In an assessment of 
the relative importance of the various terms in the equations Peterson 
(1972b) concluded that: 
(a) The horizontal momentum and turbulent energy equations are essential 
in a physically reasonable model of the flow. 
(b) The vertical momentum equation is less important. 
(c) Neglect of the vertical motion terms and the continuity equation 
does not alter the major features of the downstream profiles. 

The lack of similarity discovered by Shir has been observed 


in wind tunnel flow by Antonia and Luxton (1971, 1972). 


an ts 
Lb a is se ga 
qubeeteny att an a | 
wip ad geant) de nei pay wiht fan a 
se Meaniog. oF adh ed . 
atts 4447 br” id Ahatli ght WY cowl 
oot rmbereibh Sagas tate AOR trap: ons 30) Gheh com! et 
sawed (comp pemaih giqrattth do Tel lale sb Bea-ees opoed aa ell 
ee ates: bite G40) Pe. share ad) Ques Sater Ohare asa ott al 
caves @27 2% «.goqal G13 fs OLN! ganmet at (XO) sp@akw, oo) sighed 7 
sovel at) 3s Oth crew? alti fl qher ir) Gee Shy: rey) taney a teem 
th Si dime! ais@el Giiiinige! avew asf CA (rhe lata: ‘=eunoed 
distetvideend allt ailithed) Ligy) 16 1witeriion. an OM —odap eos 400 Show 
sb 209 CON haw evan (hemy-a dtpoes allt eer VOINE emds ea obeRe 
Udiioeiepetth oe aay) (qe ee, au eseigor 
cfd nie! agipiveibets ald haqitervyie oh (attgl ) aggxnset 
naires “h? oad wa mM utht! aj 87 Sal juehte ,@'uids tae saoegee 
ia geenteemee We ?° cdgiinavined 1G) tnt Coe jeu gales iston n 3a 
aeeret Gubigeipe Idt Gt tata) eioltay ade Ty hampered om tale odd 
val hohekewts faetet) 
Coten@een ose Grelieum grad sie!) aaut Oe, Gute snsonpaied wit Cad , 
wll aly to 16065 Alsen elewbeyaty a at 
Je 2NOgdl Baewel OF externa Sn io teuia yey elt Ct) 
AD ipe = ings was bie Genes worn fati pee? wile v1) davigs¥ {s) 
nd iaes manidanwob: yl? 0 op Ti teat: slew Git, Terke 2m BeOS 
beyneai. Sem, wait THAL yeh beanmenth vfouitimte 39 deal aT 
{98Or TERY waren bie Al ooh eA sas En “ad 


paeianee ee: 


: 


41 


4.4 An Alternative Turbulent Energy Model 
The hypotheses of Glushko (Beckwith and Bushnell, 1968) have 


been adapted by P, A, Taylor (1972) to the atmospheric change-of-terrain 


problem, The mean and turbulent quantities were related by assuming: 


ED) Aes Ky 3 
1/2 
(2) =E Le 


where Le is the mean scale of the turbulence given by 


* 
b. ='k 
(z + Zoi) 


G 
* 1/2 
and k =ka , a being the equilibrium constant stress layer value 
OLethe vatio. “oli. 
(3) K_=K 


E- M 
(4) The dissipation is represented by 


3/2 


E=E£ LP 


with the dissipation length L, given by 


Le= k, (225) 


where ee k ee 

The values of the constants ee and ke were chosen to be consistent with 
a local balance of production and dissipation in the equilibrium 
constant stress layer where t=a E, Taylor used the value 0.16 for 


a, as suggested by Peterson, The boundary conditions on E were 
a ieee an 


RE =O at 2-0 
There is no flux of turbulent energy through the ground, 

In general the results of the Glushko model were intermediate 
to those of the mixing-length and Peterson models, The depth of the 
modified region was larger than that of the Peterson model, but less 
than that of the mixing-length model, The inflection point in the 
velocity profile was present, but the change in gradient was not nearly 
so sharp as in the Peterson model, Figs. 4-5, 4-6, and 4-7 compare 


the three models (and the modified Peterson version), Details of the 


numerical methods are presented in Appendix B, Taylor also ran tests 
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with different values for a and found that larger values produced slightly 
more pronounced modifications, that is, larger stresses and deeper 
internal boundary layers, 

Panofsky and Petersen (1972) found that profiles from the 
Risg tower in Denmark were more nearly logarithmic in the modified 
region than predicted by either the Peterson or Shir models, Petersen 
and Taylor (1972), in a later analysis of the same data, indicated that 
the internal boundary layer heights predicted by the Peterson model 
were too low, The shape of the profiles appeared to be reasonably close 
to that predicted by the Glushko model, Bradley's data also suggest 


that only a weak inflection point is present, 
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, 
% a, 
ML Mixing-length model, 
4-0 \ GiL------ Glushko-type model, 
\ EP — — Peterson model, 
SAR PM—-—=Modific © Peterson model, 
\ 
VO me 
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lZ59 
(b) Rough-to-smooth flow, M = 3 
Figure 4-5. Downstream variation in friction velocity for various 


numerical models (from Taylor, 1972). 
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4.5 Other Numerical Work 

Most other numerical investigations of the internal boundary 
layer problem have employed the mixing length relation and differ only 
in the form of the boundary conditions and equations, or in the finite 
difference scheme, for example Nickerson (1967). Exchange coefficient 
theory has been employed in numerical solutions by several Russian 
workers as outlined in Panchev, Donev and Godev (1971). Smith (1967) 
used a simple K-theory approach, neglected the terms in w, and closed 
the system by developing a semi-empirical equation for the rate of 
change of K. Since K is a measure of the turbulence, it will change 
because energy is fed into or removed from the mean flow during 
acceleration or retardation, and because the turbulence diffuses itself. 
Consequently, the equation proposed was of the form 


u ak - AK OUx ae ee i 


ac9 
oO Uy az ae 


where A is an empirical constant. Although the method is an over- 
simplification, it does serve to illustrate a different line of 

thought. A general technique for the calculation of all types of 
turbulent shear flow was presented by Nee and Kovasznay (1968). 

Their formulation of the equation for the rate of development of the 
eddy viscosity included the effects of convection, diffusion, generation 
of turbulence, decay of turbulence, and straining due to acceleration 
by pressure gradients. Such an approach could, no doubt, be used to 


solve the atmospheric change of terrain problem. 
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CHAPTER V 


A MICROMETEOROLOGICAL EXPERIMENT 


DIL “Wave Experimental Site 


A Site suitable for a change-of-roughness study must meet 
several requirements: 

(1) It must be level to avoid orographically~induced vertical 
velocities, Such winds may have produced the anomalies described by 
Cameron (1970), 

(2) There must be adequate fetch across a uniform surface to 
ensure that the flow is in equilibrium up to the maximum height of 
instrumentation. The appropriate height/fetch ratio.is one of the things 
to be established in solving the change-of-terrain problem, Because it 
was virtually impossible to find any extensive uniform areas, estimates 
based on the Peterson theory were used as-a guide in the selection of a 
site, 

(3) The site must be free of large obstacles over distances 
several hundred times the greatest height of observation, 

(4) There must be a step change to a second uniform surface 
in order to reproduce the two-dimensional situation for which the 


theories were developed, 

After no little search a site which came reasonably close to 
fulfilling these requirements was located on the farm of Mr. Bill Perry 
near Chin, Alberta. A detailed map of the region is shown in Fig. 5-l. 
The section of land slopes downward to the north with a maximum relief 
of about twenty feet in a mile. The chosen site was situated at the 
edge of a mustard field on a fairly level area about midway down the 
field (east-west slope about 5 ft/mile). A uniform fetch of at least 
1/6 mile (268 m) across the mustard was available in westerly (prevailing) 
winds. In view of the preceding theories this would be adequate for an 
established wind profile at the maximum anemometer height of 8 m, but 


it might not be sufficient for equilibrium at the 4m height of the 
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fast response instruments. Immediately to the west of the mustard was a 


wheat field extending for 1/4 mile (402 m). The crops were both about 55 
cm in height, so it was expected that their roughnesses would be the 
same order of magnitude. Some of Cee ane formulae discussed in 
Section 5.5 predict exactly the same values for the two roughnesses. 
Because the mustard was in full blossom with large flowers in addition 
to its rather larger leaves and stalk, it was anticipated that its 
roughness would be somewhat greater, as the formula of Lettau (1969) 
would suggest. Consequently, a small smooth-to-rough change would be 
present. However, with a roughness~change ratio of order one, very small 
modifications would occur. Im fact, usine a 107 criterion, the flow 
characteristics of the wheat and the mustard would probably be indis- 
tinguishable. 

In the event of an easterly wind a uniform fetch of at least 
1/6 mile (268 m) was available across the adjoining fallow. The nearest 
buildings were at the farmyard over 1/2 mile away, and the nearest trees 


were over 1 mile away. 


5.2 Instrumentation 

The fifty-foot main tower was erected just in from the edge of 
the mustard field with the booms extending outwards over the mustard 
(Fig. 5-2). Normally there are six such booms placed logarithmically on 
the tower to give instrument heights of 1/2, 1, 2, 4, 8, and 16m. On 
this site the 16 m height was not employed because of inadequate fetch 
and the 1/2 m level could not be used because the mustard was some 55 cm 
Misti. . 

Mean wind speeds were measured with four RIMCO-CSIRO sensitive 
cup impulse anemometers. According to Sumner (1968) these 3 1/4-inch 
cup, low-torque anemometers have a starting speed of no more than 10 cm s- 
and are linear down to this speed to within 1 1/2 %. The accuracy is 
claimed to be +1% and to change no more than 1/24 with extensive field 


use, although starting speeds may increase typically to 19 cm s-l 
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Figure 5-2. 


Photograph showing step change from mustard to fallow and 
location of main tower. 
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Calibrations were carried out in the University of Alberta Department of 


Civil Engineering Wind Tunnel, using a hot-wire anemometer and a pitot 
tube as standards. Similar tests were also conducted in the wind tunnel 
of the Defense Research Establishment at Suffield, Alberta against a 
secondary standard, It was found that a more realistic value for the 


accuracy of these particular anemometers is+5%at speeds greater than 


1.5ms No correction was made for a possible 1% overspin in 
fluctuating winds (Hyson, 1972). 

The vertical wind was measured with a sonic anemometer mounted 
on the main tower at a height of 3.66 m. Kaimal and Businger (1963) have 
estimated the accuracy to be 154. Uncertainty in the calibration constant 
would raise this figure somewhat. The response time is about 0.1 s. 

Heat fluxes were measured with a Fluxatron (Dyer, Hicks, and 
King, 1967). The propeller has a response time of about 0.4 s and the 
temperature sensor, 0.02 s, This instrument was mounted on a separate 
mast at a height of 3.66 m and was oriented to give the best exposure for 
the wind direction on any one test day. McBean (1972) has demonstrated 
that heat fluxes determined from a propeller instrument such as this 
could be too low by as much as 25% in near-neutral conditions. 

Momentum flux was measured by the eddy correlation techni que 
described by Hicks (1969). Random errors of about +204 are expected. 
During the easterly winds of June 30 this vane-type shear-stress meter 
was mounted on a portable mast about 20 m into the fallow. During the 
westerly winds of July 1 it was mounted on the main tower. In both cases 
the height was 3.66 m. Hicks (1972) has estimated that such an instrument 
mounted at 5 m gives values about 8% too low. In view of the lower height 
used here the loss would be somewhat greater. Unfortunately the vertical 
propeller ceased to operate properly after the perspex retaining screw 
broke. 

On this same mast at heights of 1m and 2 m the two remaining 


RIMCO-CSIRO sensitive cup anemometers were placed. When downstream 
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profiles were to be taken, three of the anemometers from the main tower 


were removed and mounted on this portable mast. 

Data from the fast response instruments—the fluxatron, sonic 
anemometer, and shear stress meter— were recorded on tape by the 
multiplexing, digital conversion, minicomputer control system described 
by Honsaker, McDougall and Oracheski (1972). The mean wind was recorded 
by photographing electric impulse counters at 15 minute intervals. 

Temperature and dewpoint sensors were also mounted on the main 
tower. The former failed to function properly, although this was not 
discovered at the time. Hourly temperatures for the test days were 


obtained from the nearby Lethbridge airport. 


5.3 Experimental Procedure 


Two test days produced the required conditions. On the after- 
noon of June 30 the wind was from the northeast and complete data were 
obtained from 1300-1400 and 1600-2030 MDT. Near-neutral conditions 
prevailed on the morning of July 1 as a result of heavy overnight cloud 
cover and moderate westerly winds. From 0640-0805 MDT ten-minute wind 
profiles consisting of 5 levels up to 8 m were taken at distances of 
10, 20, and 30 m downwind from the mustard-fallow boundary. Near-neutral- 
ity was confirmed by the heat flux measurements at 0800 which gave a 
value for the Monin-Obukhov length of -770 m. Complete fast-response 
data were obtained from 1300-1600 MDT. In the afternoon near-neutral 
conditions once again prevailed (Monin-Obukhov length L=-900 m). 
Ten-minute downstream profiles were again taken for the period 1626- 
1712 MDT. Unfortunately, the wind had shifted to 320-330° , so that a 
two-dimensional situation was no longer present; effective fetches of 
22, 32, and 47 m were obtained. Uncertainty in the estimated wind 
direction could mean that these distances are in error by as much as 
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5.4 Calculation of the Wind Proffile Parameters 


In neutral conditions it is generally accepted (for example, 


Lumley and Panofsky, 1964, p. 103) that the wind profile is given by 
u ==(tt/k) In z/z Z>Z 
v oO e) 
where the symbols have the same meanings as before. 


When the surface is very rough, an additional parameter, the zero-plane 


displacement d, is often introduced so that the equation becomes 


A == C 1S 


AAD) In (z-d)/z. z>d+2. 

(for example, Munn, 1966, p. 59). 

Since its introduction by Rossby and Montgomery (1935) there has been 
considerable controversy over the exact role of the displacement d. 

Munn (1966, p. 59) stated that the modification is empirical and provides 
a better fit for experimental data. Pasquill (1962, p. 71) claimed that 
it allows for the virtually stationary layer of air trapped within the 
roughness elements. Sutton (1953, p. 239) qualified it as applicable only 
for large roughness elements that are uniform in height and distribution, 
in which case it is to be regarded as a datum level above which normal 
turbulent exchange takes place. He pointed out that if both z, and d are 
regarded as independent arbitrary constants, then the profile containing 
them cannot be derived from the first-order differential equation in the 
usual fashion. Sellers (1965, p. 150) asserted that because observed 
values are so variable and random, no real physical significance can be 
attached to the zero-plane displacement. He went on to say that in most 
cases the profile could just as well be described without the introduction 
of the parameter d. Lumley and Panofsky (1964, p. 103), Geiger (1965, 

p. 274) and others interpreted the zero-plane displacement as a shift in 
the origin of the z-axis because a vegetation cover produces an effective 
surface above the solid ground, the wind going to zero at z=z,+d, 
Lettau (1969) identified d with the corrective height increment between 


the mathematically defined zero plane and the arbitrary datum level from 
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which an observer has measured anemometer heights. Monin and Yaglom 


(1971, p. 293) derived the profile containing d formally, by treating it 
as a translation of the zero level for the vertical axis. They regarded 
it as being analogous to the displacement thickness in boundary layer 
theory. 

Given a set of wind observations at different levels under 
neutral stability, the three parameters mes Zo, and d can be evaluated 
using the method of least squares outlined by Robinson (1962). However, 
steady adiabatic conditions are infrequent; neutrality is most often a 
transient state occurring between night-time stability and day-time 
instability of the surface boundary layer over land. Usually the 
parameters must be determined from winds observed during unstable condi- 
tions. A general diabatic profile is given by Lumley and Panofsky (1964, 


pevlle)ias 
u =(u, /k)f in (z-d)/z, - vef)| 


where ¥ , a function representing the diabatic influence, is given by 
y Leff wesc ys 
6 7 
where $ is the non-dimensional shear, and § is either g/l; or z/L* 
depending on whether heat flux or temperature gradient is available. 


The Monin-Obukhov length is given by 


3 
L=-u, c, pT/kel 
whereic is the specific heat at constant pressure, O is the density of 
air, T is the mean air temperature, k is von Karman‘s constant, H is the 
heat flux, and g is the acceleration due to gravity. The related length 


L' was introduced by Panofsky (1963) as 


t — 8 i) 
iF ie 2 T / ok 26 


where © is potential temperature. Therefore LI=K/Ky L , where 


Ky is the turbulent transfer coefficient for heat and Ky is the turbulent 
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transfer coefficient for momentum. Both parameters are related to the 


Richardson number, eas 


such toate, Lik 1 p : 
An equation which interpolates between neutral stability and 
free convection has been advanced by a number of authors (Lumley and 


Panofsky, 1964, p. 110), namely, 
3 
pi -—ye oo =| 


where ¥ is an empirical constant. 


Substitution of b=(z/L') /Ri yields 


Zz Ri 


— OO eee 


2 t 
a Gare): 


and elimination of z/L' gives 
Sane 


$ = (1- ¥ Ri) 

This is the so-called KEYPS formula. The value of ¥ is usually 
determined by numerical integration. An alternative formulation by 
Webb (Panofsky, 1963) involves a two-part explicit function for W . 

Paulson (1970) provided a mathematical representation for ¥ 
by integrating the analytical expression for the non-dimensional shear 
obtained under the Businger-Dyer hypothesis that z/L=Ri. This hypothesis 
was based on the Kerang observations (Swinbank, 1964) and has also been 
suggested by Pandolfo (1964). An immediate consequence is that 
Ky/Ky =1/0 . Hence the ratio is one in neutral conditions and increases 
with decreasing stability, in agreement with a suggestion by Priestley 
and Swinbank (1947). It also helps to account for the disparity of 


reported values for Ky/Ky (lumley and Panofsky, ‘1964, p. 105). 


‘ of — V4 
With $ = Ce | 


the diabatic influence function becomes 


oe | ' 
ee min 


i , 7 7 > i > J a | : o 
tow gh OD IRS tetesoen, pete’ Lee ada 
tee ge duad) poastn vo vata 6 ee a in oe mnie s, 


‘» fp ty 


Semcon Look hye aml Bey 2 
oehekg . AA ¢'ys)=4, Gan 


4a 22 Nie 
OT sti - ; 


Simitihot = ta Gaped at » nied CON bel her-oe ri: 

qf oofseiamant aviteta te Gk: «eebaereeee fesce tire vd beat 

. © 2) wehonil 101 (gk daapoh © eeuteomh (ag! vpledone®) #2 
=) plvitt=.w) gs leshtoontigh — i hphigneg (AEs) enelant | 

eet) jiu bea. aly Se. eoheaattee Leskonehane, oe gntang ; 

disudsgged ¢i€f (Ne SW soaks oleate 29/0 reqpehaiatt ents pare 

mot velo aud hes 1A) dimdnred) secant gaeeet oar aa 

jet) ¢) ose Goel @enhip> BA bgbt) oifelun’t qd be: } Bs 

ACanneee) (we arelelhnos betauay ub ano ut Ohiagd aga saneh a Wt = 

eflaate) q@ dalewugeva 6 late dimmrrge wh gene tee ptianetoah 2 = 

Hapa et oie Seed acumen ie 

«QOL 9g 08Ok pean ireatt Bem. <a : 


ee aad 
Pee ve nee tS ridgias 
a > aes 7; 


56 
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/ 
LENG + = - 7 
Cea Ome (Pe + bot 4] a 
i 
where r=(1l- ¥z/L) - 
This implies that for -z/L large u varies asymptotically as hs in 
-1/3 

contrast to KEYPS where u%z f . From analysis of the Kerang data 


Paulson found the value for ¥ to be 16. 

Using this representation, the wind profile parameters were 
computed by the method of Stearns (1970). The roughness and displacement 
height are determined so that the sum of error squares between the 
observed and calculated wind speeds is a minimum. The friction velocity 
is treated implicitly as a function of d and z,. ‘If &; is the error 


between the measured wind speed and the theoretical wind speed at height 
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Successive approximations to Z and d were obtained using the secant 


method, With an initial estimate fo. d and u, two guesses were made for 


Zz, and approximations were obtained from 
z(n) GZCne) ~ Zens!) © 2En) 
z,(n+2) = ? 
GA nti je, Ge Cn) 


until the difference between successive approximations was less than 
0.01 ecm. This value for Zo together with d(1) and a corrected uy, were 
then used to compute GD(1). A second guess d(2) was made, and the 
procedure repeated to obtain a value for GD(2). Successive estimates 


for d were then obtained from 
et Gy) en chen = (wei) Cte) 


d(n+a2) = aa Z ae 


until the difference was less than 0.1 cm. 

The wind-profile parameters for the mustard were computed from 
twenty hourly-wind averages taken at four levels on July 1. For the 
fallow, d was assumed to be zero and the parameters were computed from 
twenty hourly-wind averages recorded at the two levels on the portable 
mast on June 30. The data are displayed in Appendix A, 

The value of uy under neutral conditions can be related to the 
mean wind at some level by means of a drag coefficient, namely, 

Cq =U, /U} , where Ca is the drag coefficient, and Uy is the wind -at | mm, 
Bradley (1972) has shown that even in unstable conditions, such a relation- 
ship still holds for low levels, in keeping with the commonly-held view 
that neutral conditions always prevail at sufficiently low heights. 

The results of the computations appear in Table 5-1. The 
roughness ratio for the site is thus m=12:1 and the roughness-change 


parameter is M=2.5. 
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Table 5-1. Results based on twenty hourly-profiles, 


Surface Roughness Displacement Drag Coefficient 
vs SiGe d fs) Sieh fe Sisfis 
fe) d 
(cm) (cm) (cm) (cm) 

mustard 12 eC 15 led 0.204 0.008 

fallow Oe Oe? 0 --- 0.091 0.004 


Table 5-2, Some values for the zero-plane displacement of mustard 
55 cm high and similar crops, 


Value Circumstances Source 
(cm) 
jis) mustard 55 cm high observed at Chin, Alberta 
55 d=h., from models in a 
wind® tunnel Plate and Quraishi (1965) 
27-55 Jal b/s cont slv=aiet Monin and Yaglom (1971) 
10 grass 50 cm high Taylor (1962) 
Pg grass 60-70 cm high Deacon (Geiger, 1966, p. 274) 
30 grass 60-70 cm high Calders (subton,. 19554 9.0240) 


55 d=h), for field crops Paeshke (ibid., p. 239) 
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5.5 Appraisal of the Computed Parameters 


Values for the displacement height are not plentiful in the 
literature, but some empirical formulae do exist. Sutton (1953, p. 239) 
quoted Paeshke's result for dense crops that d is equal to the average 
of the measured heights of the roughness elements. Plate and Quraishi 
(1965) found for model crops in a wind tunnel that to a good approximation 
d=h, where h, is the height of the crop. Both experiments also indicated 
that the velocity profile is logarithmic only for 272005 s0-calied 
‘canopy flow' governed by the nature of the individual roughness elements 
occurs below this height. Taylor (1962) chose a value Of m= 5 bh for 
longegrass, 50 cm high. Plate (1971, p. 29) stated that in general d will 
differs from he when the) density of the roughness elements is sparse. In 
the limit the flow may not be fully rough and special investigation is 
necessary. Geiger (1966, p. 274) asserted that d is not equal to the 
average height of the vegetation surface and that only an examination of 
the wind profile can determine the value of d which will best satisfy the 
equation. Geiger also quoted some results of Deacon where a displacement 
of 25 cm was used for long grass 60-70 cm high, that is, d=2/5 ho. 

Monin and Yaglom (1971, p. 294) took d to be zero for low vegetation 
and suggested the range h/2 adh) for high= vegetation. 

Geiger attributed this wide range of values to the high degree 
of pliability of the plant stalks. Rauner (Munn, 1966, p. 161) found 
that as the wind speed increased over a forest the value of d decreased 
from h /2 to h/4. 

The value obtained for the mustard was 15 cm. Even though the 


' a rather large stan- 


wind speed at 1 m was quite steady at 2-2 1/2 ms" 
dard deviation of 4 cm was found. The ratio d/h, for the mustard was 

3/11 or 0.28. This was lower that most estimates, but higher than that 
of Taylor (1962). A comparison of values computed from empirical form- 


ulae and those reported by various workers for similar crops is given 


in Table 5-2. 
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No micrometeorological studies appear to have been done over 
mustard; there are no published values for the roughness, There are, 
however, once again a number of empirical formulae, Plate and Quraishi 
(1965) found that for crop-like elements in a wind tunnel Zed We, 
in agreement with the field work of Paeshke (Plate, 1971, p. 27). 
According to Sellers (1965, p, 150) Tanner and Pelton, Kung, and others 
have used relationships of the form 


log z= log a+b log he 
thatecs: 2. =a ne 
O c 


where a and b are empirical constants, Values of b between 0,99 and 
1,42 have been employed and values of log a ranging from -1.4 to -0, 88, 
Plate (1971, p. 29) objected to an equation of this form because it is 
not dimensionally consistent, Monin and Yaglom (1971, p. 294) 
recommended that h /10¢zi< ho/5; noting that the constant of 
proportionality is considerably greater than the value 1/30 for 
Nikuradse's sand roughness, They attributed the variability of estimates 
for both d and Zo to the fact that both parameters depend on fairly fine 
details of surface structure, Natural vegetation, even of one type, 

is not likely to have identical structure in different places around 

the world, 

Lettau (1969) held that the use of only the height in a 
roughness formula is too restrictive and produces oversimplified 
results, The shape and spacing of the roughness elements must also 
be considered, On the basis of the bushel basket experiments he 


proposed the relation 
Dah Sf st fe 
fe) Cc 


where s is the silhouette area of the average obstacle, that is, the 
cross-sectional area 'seen' by the wind in an approach to a typical 
roughness element; S is the specific area or lot area, that is, ,the 
average area on the earth's surface occupied by a single roughness 


element, If N elements are located on a site of area A, then SiN: 
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The factor 0.5 corresponds to the average drag coefficient of the char- 


acteristic element. Lettau considered estimates from this formula to be 
within *25%. Its validity is restricted to situations were s<S$; other- 
wise serious overestimation occurs, 

Table 5-3 gives a comparison of the value obtained for the 
mustard (12 cm) with values for similar crops reported in the literature, 
and with values computed from the various formulae. Table 5-4 shows 


Values for the fallow. 


5.6 The Downstream Profiles 

The profiles observed during the morning and afternoon runs of 
July 1 are shown in Fig. 5-3. The velocities were scaled with respect 
to the wind at the highest observation level of 8 m, which is outside the 
modified region for all fetches according to the Elliot theory. The 
uppermost points of each profile lie on the line characteristic of the 
upstream flow, except the 4m point at 10 m fetch which seems to be 
somewhat in error. The acceleration of the flow at increasing distances 
appears as a shift to higher speeds in the lower portion of the profiles. 
The number of points lying on the upstream profile decreases with distance. 
Qualitatively, these features are similar to those observed by Bradley 
(1968). 

The height of the internal boundary layer cannot be well 
defined by so few points, so direct comparison with predictions from 
the various theories has not been attempted. The surface friction 
velocities scaled with respect to the upstream value were computed from 
the 1/2 m wind speed by means of a drag coefficient similar to that 
described in Section 5.5. Assuming that the 1/2 m wind is in the portion 
of the profile characteristic of the new surface 


u(l/2 m) = (ay./k) 1n750/1 


Uso = Cao u(1/2 m) 
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Table 5-3, Some values for the roughness of mustard 55 cm high and 


similar crops, 


Value Circumstances 
(cm) 
12 mustard 55 cm high 
8 an tla) fh j2rom models 
in a wind tunnel 
Pas wheat 60 cm high 
11-15 grass.6Q-70 cm high 
cee z shit /12 
oe We 
5,9-L1 ho /lO<2 2h 1/5 
c Gr < 
6-9 long grass 60-70 cm high 
os thick erases wp to 50 cm 
4-5 high grass, wheat 
4-9 long grass 60 cm high 
14 fully grown root crops 
4-12 tall erass,. erain 
2 grass 60-70 cm high 
14 roughness element formula 


(s~20 cm2, S~40 cm?) 


Table 5-4, Some values for the roughness of 


Value Circumstances 

(cm) 
LO fallow with some trash 
Pees | fallow field 
0,6-2,0 prairie grass, countryside 
Ooms 5 mown grass 1,5-4,5 cm high 
Os7l=2. 3 downland, thick grass up 


to LOMchighieh 
5 short grass 
Ae grassy surface 
i flat country 


Source 


observed at Chin, Alberta 


Plate and Ouraish? (1965) 
Penman and Long (Sellers, 
1965, p. 151) 
Deacon (ibid, ) 


Lettau (1967) 
Monin and Yaglom (1971) 


Priestiy (1959.0, 22) 

Sutton (1953= p. 233) 

Hess (L950 3 mre 277) 

Pasquill (1962. p, 72) 

ibid, 

Geiger (1966, p, 275) 

Calder (Sutton, 1953, p. 240) 


Lettau (1969) 


eel iow? slave Sima leis Siwieei ia, 


Source 


observed at Chin, Alberta 
Hesse (19590 >, 277) 
Geiger (1966, p. 275) 
Priestly (1959, p22 


burton (1953595. 22.2) 
Calder (ibid., p. 240) 
Plate (1973-45, 27) 
iMaplel 
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where Cay — k/ln 50/1 


The wind at 8 m is characteristic of the upstream flow, 


u(8 m) = (ux, /k) in (800/12) 


assuming negligible streamline displacement. 
Uy = cq u(8 m) 

where c4,= k/1In (800/12) 
Hence Uo /Uy1 = (cgo/cg,) bua/2 m)/u(8 m)]. 
Figure 5-4 compares values of the friction velocity computed in this way 
with three analytical solutions, Figure 5-5 does the same for three 
numerical solutions. Details of the computational schemes are in 
Appendix B, The initial overshoot was greater and the recovery was 
rather more rapid than predicted by any of the theories. At larger 
fetches the Glushko model seems to be in closest agreement with the data. 
Extrapolation of the curve through the observed points suggests that the 
Peterson model might come into closer accord at greater downwind distances, 

To facilitate direct comparison of the observed velocity profiles 
with the results of the various theories, a transformation was employed 


to scale the speeds with respect to the upstream friction velocity, 
u/u,, = l/u(8 m)] [us nin! 


The first factor is the observed value and u(8 m) /U, 5 can be calculated 
because the 8 m wind is characteristic of the upstream flow, that is, 


u(8 m)/u.. = Clk) ine (800712) 


Comparisons of the observed profiles with those predicted by the various 
theories are shown in Figs. 5-6 through 5-17. All of the theories 

Aver eecinate the speeds at the shorter fetches. One reason for this 
lies in the zero-plane displacement associated with the mustard. No 
theory takes into account the possible effects of such a displacement 

on the flow modification. At 10 m downwind from the discontinuity the 


flow had not yet descended fully onto the fallow; this is indicated by 


7 
(work meoningy a7 Ye otited eetan meld ‘a 
paee nl Opeth te te” 
saan taht vot tamete 7, 
to Wings = we : 
(Li\OGS> afl ty) a7 
(ime Bum SVED a, s\n dl ee yeweee 
yaw @hal ) lade Vitwuloy ae) isind Ses, s4 LiKe germqnos Dok sTUQhT 
ay ac) are of) GaN Ct oayght anad tttas pane ‘sods at an 
14 @fnnaze Innoi®esuenas att 4h Clit veawlee ‘teska ra 
“Oy Wier oy’ lida {petaen ininigtev) Selsaiat ae? A wth . 
gvel SA, (ant veRRt al) 30 var ¥4a'i ie hete ats bhyes, st6m sat 
emad ofa citie ape jeady teas to nt ad <9 Ompde Eetess oofifeel® ol? 2 fon 
ois Init etpoygite eliiay tmv Topo ail fgvaels avis ee Ou eobent 
' balwinwwh 4osnete 2p baosar TAdlo ese ate 2 by hei Ls sor 
esliterny & ) Lavtenbhs ot) 26 noaktaes Sag7th ‘nthiians af 
orynrqe chy hokmimetan 7 » au ep awiltey sef ain tolueen nas bi a 
MPIoisy FOP 41? aaa O17? td desgege ite whoge oid =. 
hy a ir Padi Te Med = yw a 
bate {por (mes thier Oi tem oleg Ged niet — rt yonant 1” 2 odt 
cob Daal? ail on gl sit Ye abies! soiea ‘Arta 4 Arh ne ad? 
CEIYOD8) al ad gee ini 
serra Gd Yt bosulta@asy Sno r ot ess bern a 
esleneds iadiae iA de iguored +e swale 
oh) wi nGeaes Sut <aaiaan, ea f on 
ol olrmebetis sind aay i arls eaal mening sae who 
tren tutes a wie \a lendiite onda 
sé eniambanasat = oe Aves eager A eT ¢ 
wv bears Wi ak de a ee a wtiui teh 


enn reer). 


Of 


0.60 


Ux?2/Ux1 


0.50 


0.40 


0.00 


Figure 54. 
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10.00 20.00 30.00 4u0.00 50.00 
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PM modified Peterson 


GL Glushko 
ML mixing-length 
oO observed 


Downwind variation of surface friction velocity inferred 
from the wind at 1/2 m compared with three numerical 
solutions. 
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the fact that the’ upper four” points of the profile fit best to a tine 

that is based on a non-zero displacement length. A 20 ma line based on 
the roughness length alone fit the upper three points of the observed 
profiles, The flow descent seemed to be complete. However, at heights 
greater than 2 m there is little difference in the lines. Another 
possible reason for the discrepancy is the proximity to the discontinuity. 
Elliot's theory specifically applies at distances greater than 10° Z 59? 
and all other models except that of Shir exclude the region where local 
pressure effects may be felt, 

There did not appear to be an inflection point in the 
observed profiles, However, no firm conclusion could be reached on the 
basis of so few data points; a larger number would be necessary to detect 
a change in curvature, The lower portions of the profiles were more 
nearly logarithmic than the Peterson model predicted. This is in agree- 
ment with the findings of Petersen and Taylor (1972). Both the Elliot 
and Glushko solutions provide a reasonable fit to the data, In view of 
the surface shear stress calculations the Glushko model seems to be the 


best overall description of the flow modification, 


5.7 Some Turbulence Statistics 

Existing data for the moments of the turbulent fluctuations 
are inexact, strongly scattered, incomplete, and do not agree well with 
each other, The explanation is evidently the considerable error in 
measurements of turbulent quantities that have quite large variability 
both in the vertical and the horizontal, even above a relatively 
homogeneous surface (Monin and Yaglom, 1971, p, 518), The fast response 
data obtained at Chin made possible the calculation of the standard 
deviations of the vertical and longitudinal wind components, The vertical 
wind was measured by the sonic anemometer, The wind speed was measured 
by the vane-mounted horizontal propeller of the shear stress meter, and 
the standard deviation of the longitudinal component can be obtained 
from this if it is assumed that the statistics of the two variables are 


interchangeable (Lumley and Panofsky, 1964, p, 154). 
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The characteristics of the vertical component near the ground 
are reasonably well understood; the theory is relatively simple 
because the proximity of the surface prevents low-frequency oscillations 
from forming, Furthermore, the normal roughness length serves to rep- 
resent the terrain effect on the intensity of vertical turbulence 
(Lumley and Panofsky, 1964, p, 131), There is no doubt that large-scale 
features are unimportant (Monin and Yaglom, 1971, p. 518). In neutral 
conditions the standard deviation of the vertical velocity is proportional 
tothe friction veloeLty,ijthat its; 


C. wy 


Ww 


The value of the constant A, however, is a subject of much controversy, 
Estimates have ranged from 0,7 toi1.4, The mean value of all estimates 
from atmospheric observations lies in the range 1,2 to 1,3 (Monin and 
Yaglom, 1971, p., 520). A short summary of estimates is included in 
Table 5-5, Monin-Obukhov similarity theory predicts 

By /u, =F (Z/L) 


where Fy is a universal function of the stability parameter z/L, The 
nature of the variation with stability is likewise a subject of debate, 
Recent work suggests that the ratio is not dependent on stability except 
at large instabilities where a slow increase is noted (Bowne and Ball, 
1970. Monin and*Yaglom, 1971; p59). 

The standard deviation of the longitudinal component in 


neutral conditions is also proportional to the friction velocity, 


Oy = Buy, 


Estimates of the constant B have ranged from 2,1 to 3,0, Recent 
atmospheric data suggest a value 2,5 (Bowne and Ball, 1970), while 
laboratory work suggests the value 2,3 (Monin and Yaglom, 1971, p. 280)... 
A short summary of estimates in included in Table 5-6, There is evidence 
that the value varies with terrain (Lumley and Panofsky, 1964, p. 155), 
Large-scale factors may well be important in addition to roughness, 
Again similarity theory predicts 

Sa /u, =F, (z/L) 
Variations with stability tend to be small, with perhaps a slight 


increase with increasing instability (Monin and Yaglom, 1971, p. 518), 
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Table 5-5, Summary of values for ratio Ow fu. 


Value Source 

E25 Panofsky and McCc :mick (1960) 

oy Pasquill, (1972) 

Oo7 Gurvich (Lumley and Panofsky, 1964, p. 135) 

0, 87 Perepelkina (ibid.) 

Eee Klug (1965) 

Le? Panofsky and Prasad (1965) 

ee Mordukhovich and Tsvang (Pasquill, 1972) 

133 Businger et al (1967) 

1e29 Busch and Panofsky (1968) 

pie b) Haugen et al (1971) 

1.4 McBean (1971) 

Cag Laufer, Klebanoff, Comte-Bellot, Coanti 
(Monin and Yaglom, 1971, p. 280) 

0.9 Cermak, Sandborn, Chuang 


(Monin and Yaglom, 1971, p. 520) 


Table 5-6, Summary of values for the ratio 6y /u,. 


Value Site Source 

Beas various Davenport (1961) 

29 O'Neill Lumley and Panofsky (1964, p. 155) 

ie Australia saw iie | 

Parca Brookhaven ibid, 

Lo Tsimlyansk Monin (ibid, ) 

ade) pipe flow Dodds 

on) O'Neill Klug (1965) 

Zea laboratory Laufer, Klebanoff, Comte-Bellot, Coanti 


(Monin and Yaglom, 1971, p. 280) 
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Values of ©, and © were computed using a 15 minute averaging 
period to coincide with the nominal 15 minute recording interval of the 
mean wind speeds, Four such values were then averaged to produce an 
hourly figure. This procedure was followed by Haugen, Kaimal, and 


Bradley (1971). The corresponding hourly values for u, were then used 


to calculate the ratios SO ./u, and su/u,. Twelve ratios on June 30 gave 
Cwl/u,=0.94 S eyes Oo 
Cal ate AP Sieg Out 


The stability for the period of these measurements ranged from z/L= -0,39., 
t0 z/L= -0,080. Ten ratios on July 1 eave 
Cw/u,= 0.96 sod. ,O0,.18 


Suit 2 2 sod.4 0,2 


with stability ranging from z/L= -0.075 to -0,003, The ratio oye ibe 
exhibited no significant difference for these two broad stability 
categories, If the indications are correct the marked difference. between 
thetvaliues for Cu lu, must be attributed .to the effects of terrain, 
Although no obvious large-scale differences were observed, a topographical 
map disclosed that the northeast flow of June 30 was subject to different 
large-scale influences from the westerly flow of July 1, The average 
value for theeratio Ow /uy, is 0.95, which is in excellent agreement with 
the value for laboratory flows, The average value for the ratio 6, /u,, 

is 2,6, in fair agreement with other recent atmospheric values, 

Since E31 /2)Cou 6s Pos ) eaevalue (a0n thewcoustantyot 
proportionality between shear stress and turbulent energy can be obtained, 
provided the ratio 6,/u,, is known, No measurements of the standard 
deviation of the lateral component were made, As with the other statistics 
the dependence of Oy Ju, on stability and/or terrain is an open issue, 

In neutral conditions Monin and Yaglom (1971, p. 518) gave the value 

1.7, the same as for laboratory flows, Lumley and Panofsky (1964, p, 

146) gave values ranging from 1,5 to 2,6, Takeuchi (Bowne and Ball, 

1970) reported the value 4,0, while Elderkin (Bowne and Ball, 1970) 
obtained the value.2.0, Assuming PED Pome 8) Ug) together jwith the 
average observed values of ‘w= 0.95 u, and o,4 = 2,6 u, the result 
the constant lying between the value 0.16 used by Peterson 


is E=0,18 Ga 


cies 


and .the value 0.22 used by Shir. 
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CHAPTER VIL 
MORE TOWARDS REALITY 


Study of the change-of-roughness problem originated with the 
desire to understand more fully the effects of real terrain on the 
atmosphere, Substantial effort has been expended and some progress 
has been made, In developing the theory certain assumptions were 
necessary to render the problem tractable, The most common assumptions 
were: 

(1) neutral stability 
(2) two-dimensional, singlé, abrupt change of roughness 

(3) change in roughness alone, no change in other variables, 
These assumptions effectively restricted the theory to somewhat 
unrealistic situations, For a model to be an adequate representation 
of the interaction between the Earth's surface and the atmospheric 
boundary layer, other, albeit complicating, factors must also be 


considered, 


6.1 Non-Neutral Stability 

Because the theory is restricted to adiabatic conditions, 
experimental testing is quite tedious, Neutral stability, as noted in 
the preceding chapter, is rather infrequent, Awaiting the time when 
such conditions do present themselves is an onerous task, Onshore flow 
at a coast is an approximation that has been the source of much of the 
data at present. For the surface boundary layer over continental land 
areas the extension of the theory to include all stabilities is 
imperative, 

Elliot (1958) reckoned tentatively that instability leads to 
a small increase in the height of the internal boundary layer, and that 
stability leads to a small decrease, The observations of Echols and 
Wagner (1972) supported this idea, and suggested further that the 
internal boundary layer height is dependent on wind speed, In near- 
neutral conditions the height decreased with increased wind speed as 


would be expected from classical theory of boundary layer growth over 
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a flat plate (Schlichting, 1968, p, 599). In stable night-time conditions 
the opposite occurred: the height increased with increased wind speed. 
The reason, they surmised, was a tendency toward laminar flow, which 

has lower boundary-layer heights according to the classical flat-plate 
theory (Schlichting, 1968, p. 130). Increased wind speed increases 
mechanically-induced turbulence and destroys the laminar character of 

the flow. 

Blackadar, Glass and Panofsky (1967) developed a perturbation- 
type theory employing the log-linear wind profile for non-adiabatic 
lapse rates, Integration of the equations of motion across the 
modified region, with the constraints of continuity on the stress, 
wind speed, and wind shear at the top of the layer, produced an 
expression for oh tel VAL) as a function of the downwind distance 
from the discontinuity. The results were compared with the theory of 
Miyake as modified by Panofsky (1967) to include diabatic effects. 

Both theories predicted an increase in the internal boundary layer 
height with increasing instability, but the latter indicated a larger 
effect, Their experimental data was too scattered to discriminate 
between the two theories, 

Several of the other models also lend themselves readily to 
extension into non-adiabatic conditions, The analysis of R, J. Taylor 
(1962) can be made applicable simply by using the streamfunction 
appropriate to the diabatic wind profile, The treatment of Dyer (1962) 
could be modified by making the exponents in the power laws functions 
of stability. Once the variation of the ratio t/E with stability is 
known, the Peterson and Glushko models can also be adapted. 

Pasquil46(1972) noticedathatsthe cunverion the ginteriace 
height is similar in shape to the curve for the increase with distance 
of the mean vertical displacement of passive particles released from a 
crosswind line on the ground, He drew a simple analogy between mass 
and momentum transfer, the 'source' of passive particles being identified 
with the 'sink' of momentum, By using the gradient transfer approach 
for vertical diffusion and adopting realistic profiles of eddy diffusiv- 
ity and wind speed, he was able to extend the internal boundary layer 


calculations to conditions of non-neutral thermal stratification. As 
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before, the height was found to be greatest in superadiabatic and least 


in stable lapse rates. 


6,2 Simultaneous Changes in Other Elements 


A further complication arives when there are, in addition to 
a change of roughness, simultaneous changes in other meteorological 
elements, Although the mustard-fallow transition was treated simply 
as a change in roughness, in reality, there were also discontinuities 
in surface temperature, heat flux, and humidity. The assessment of 
Dyer (1962) was intended to provide working estimates for the mod- 
ification of these profiles and fluxes, It does not, however, allow 
for any interactions between the diffusing agents. Some progress along 
these lines has been made by P, A. Taylor (1970, 1971). 

The transformation of other meteorological variables besides 
the wind profile has been studied by Nadejdina and Novikova (Panchev 
et al, 1971), using the energy balance equation and the Kolmogorov 
expressions for the exchange coefficient and dissipation rate. An 
approximate solution for the transformation of the wind profile alone 
as a result of a simultaneous change of roughness and temperature 
stratification has been obtained by Dmitriev and Sokolova (Panchev et 
al, 1971), using the wind profile-exchange coefficient model of 
Laikhtman (1964). The solution revealed that the influence of the 
roughness change is less than the influence of the change in intensity 
of turbulent exchange due to the different stratification in upstream 
and downstream flows. This enabled Zaitsev (Panchev et al, 1971) to 
simplify the equations to the case of a change in stability alone, 


for which he was able to obtain an exact solution in Bessel functions, 
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6,3 Fatchy Roughness 

Yet another complication is the finite lateral extent of 
different surfaces in real terrain, The theory needs to be extended 
to include the three-dimensional case, In addition, real terrain is 
seldom simple enough to be represented by one or even several abrupt 
changes from one uniform surface to another, Internal boundary layer 
calculations become futile. Smith (1967) represented complex terrain 
by allowing the drag coefficient to vary about a mean value with a 
specified standard deviation. The results for a periodic variation 
showed that the height to which flow modifications penetrate depends 
essentially on the period of oscillation, and to a smaller degree on 
the amplitude, An observer flying low over an area made up of fields 
of different roughness would record the oscillations, but the higher 
he flew, the smaller would be the variation, The effects are damped 
with altitude, 

Pasquill (1972) applied his 'particle source-momentum sink' 
analogy. A particular elementary area in isolation, he reasoned, is 
a 'sink' which will produce a momentum deficit at any given height that 
will rise to a maximum at a certain distance downwind and then fall off 
continuously as the distance is increased further, The functional form, 
he noted, is the same as the ground-level distribution from an elevated 
source, Using this it was possible to define the upwind area of ground 


' dominate the momentum deficit at any 


within which the unit ‘sinks 
particular height. The crosswind and alongwind dimensions of these 
effective areas were derived from estimates of the vertical spread. 

The dimensions are fairly sensitive functions of height and are further 
affected by stability and roughness in that order, These areas are not, 
of course, the only roughness elements affecting the flow properties at 

a point. They are, however, the dominant ones, Any outstanding bigger 
feature outside the area may also exert an important effect, Nor is 

this method rigorous, In order to obtain quantitative estimates, 

Pasquill was forced to assume that the diffusivity field is horizontally 
uniform, which it patently is not. It must reflect the spatial variations 


of the very patchiness in question. Nevertheless, he felt that the 


results would be useful as a first estimate in evaluating single point 


measurements, 
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Accordingly, measurements made at some position will be 
representative of the underlying patchy terrain only if the patchiness 
is on a scale small compared with the effective source area, As a 
concrete example, in flat rural England the scale within which most 
of the major roughness elements are ‘epeated is about 100 m. Pasquill's 
calculations indicated that measurements at a height of 50 m would 
probably be representative, except perhaps in unstable conditions, 

In areas containing cities with tall buildings the scale of the varia- 
tion would be some kilometers, and the height required for representative 
measurements, some hundreds of meters. Attention is again focussed on 
the use of low-flying aircraft as an approach to direct spatial averag- 
ing. Fiedler and Panofsky (1972) deduced that the effective roughness 

of the Pennsylvania plains is about 40 cm, considerably larger than the 

1 cm figure usually ascribed to grassland. The importance of individual 
obstacles ranging from bushes and hedges to woods and buildings is 
thereby confirmed, 

The roughness formula of Lettau (1969) permits the objective 
determination of the dominant contributor to the effective roughness 
among such a variety of roughness elements. The significance of the 
natural topography of the landscape is still largely unknown, ie 
Lettau's formula can be applied on the continental scale, the effect on 
the aerodynamic roughness can be assessed, Another possibility is 
spectral analysis of the total variance given by a topographical 
profile curve. Such a technique would numerically express the relative 
importance of various topographical features—knolls, ridges, mountains, 
etc, -— qualified, of course, by their geographical frequency. It should 
therefore be possible to extend the roughness change theory into the 
realm of the planetary boundary layer, The effect of the whole spectrum 
of roughness on the geostrophic drag coefficient could have important 
repercussions on the atmosphere as a whole, Some work in this area has 


been done by P. A. Taylor (1969b)., 
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CHAPTER VII 


CONCLUSION 


TEL Summary 


The zero-plane displacement for mustard 55 cm high was 
found to be 15 cm and to exhibit rather large variations (standard 
deviation 4 cm), This value was lower than might be expected from 
most empirical formulae, The roughness of the mustard was determined 
to be 12 cm, in close agreement with several estimates, most notably 
the roughness element description formula of Lettau (1969) and the 
power law quoted by Sellers (1965, p. 150). The fallow was found to 
have a roughness of 1 cm. 

Analysis of the downwind profiles revealed several things. 
The surface friction velocity inferred from the wind at 1/2 m overshot 
initially and returned slowly to equilibrium, The wind profile in the 
lower part of the internal boundary layer was more nearly logarithmic 
than predicted by the Peterson model, In other words, the height of 
the internal boundary layer was somewhat larger than predicted. No 
inflection point was present in the wind profile, Although the Elliot 
theory was in fair agreement with the data, overall, the Glushko model 
appeared to be the best description of the flow modification, The 
presence of a zero-plane displacement complicated the issue, especially 
as the theories do not take it into account, 

A plausible value of 0,18 for the constant of proportionality 
between shear stress and turbulent energy was obtained using the 
observed average values of Sw /u,=0.95, SG/uy,=2.6, and the assumed 
value 6, /u,=2.0, The two test days represented two broad stability 
categories, No variation of oy, /u, with stability was found, However, 
there were significant differences in the values of Cu alas This 
could have been the effect of stability, or, because the wind blew 
from different directions on the two test days, this could be attributed 


to the effect of large-scale properties of the terrain, 
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7.2 Recommendations 

In weighing current knowledge a few general facts may be 
Stated in addition to those directly disclosed by the experiment. 

The height of the internal boundary layer, regardless of its definition, 
grows with the 4/5 power of distance from the discontinuity. For 
representative wind profiles a height/fetch ratio of 1/20 seems to be 
adequate, For the turbulence statistics to be representative it 
appears that a height/fetch ratio of 1/100 to 1/200 is necessary, 
depending, perhaps, on the sense of the change. The turbulence adjusts 
rather more slowly than the mean wind, 

Even for the simple two-dimensional change-of-roughness 
several practical difficulties remain. The solutions are fairly 
sensitive functions of the roughness ratio, In most applications the 
downstream roughness is not easy to obtain. At Chin, reversal of the 
wind direction made it possible to calculate the fallow roughness from 
observations at the same site as for the mustard. In many instances, 
however, the instrumentation may well have to be moved in order to 
obtain representative profiles of the second surface from which to 
derive the downwind roughness, The roughness formula of Lettau (1969) 
furnishes a means of estimating the second roughness, though with poor 
accuracy. An improved method of roughness estimation based on roughness 
element description would be most desirable, The zero-plane displace- 
ment is another issue. None of the theories have included this parameter 
despite its obvious importance when dealing with forested or urban 
areas, Existing models will have to be revised in some way to deal with 
this aspect of the flow situation, New models ought to contain it from 
the outset, A further difficulty is the possible dependence of both 
the roughness and displacement lengths on the wind speed, 

No definitive final solution is yet available, even for the 
simple change-of-roughness problem, One obstacle has been the paucity 
of high quality data. Future experiments should employ as many 
anemometer heights as possible, particularly in the controversial 
blending region, to detect any change of curvature, Some effort should 
also be made to measure vertical profiles of shear stress as well as 


the downstream variation of the surface value, Other research needs 


at vaw ogee! Tegaaag we) v oplulteinad suet slgien 8 a 7 
ramizegee alle yd bea-als ah pianet hts soit eo ons adel  besese 
oowsinten «il So-waplvenger coat yralamod. tuersrat oly Youadpted adf 
a vy luepueoadth fe eos! saneedh de sailg chi mata sts Motg 
‘il vane OTE 96 stir Worel\otgdoté 691)? arg tnt besa or 
} evésmineesiqas @€ ws baidedi ete a.nei cove! ada yo4% ser aupet 7 
wanaswoen 02 HOC? wm WUI\E to aeon? Aasat) 7a4q ied « ath one 
Mba wonelew’rn) ent w“arnds oft) 1D 40a GRI Gu arr he 
ee 
pitowG~ te saas (snvlenmsdibioy® ofgaia a oH, weve 7 
yvialel Gia .dimtjolua as jit 4 poGei? a0 Iealpsuq iareves 
ily evo Lbasl Faw odin art ick eh Ue: BR. Fo ' whsena) # ttivone 

lemianr wil fae o) ween $c2 0) eueqmquar ansrjacwob 
moe) ele indce het we juited i? si@enene |) oban. «© topett, bane 7 

TL? ee ee sey wmed fh) Ls smusecrsyge | 
a ! hile ojomiGns ott ae 
pelt sadn Hoprcige bale | sry oetee meoriyet eae 
(HCl) wogde! pe shite eemniies : o> Siew hnrvewsd etl sv iva 
1Ode Tt Agios . deem | meee ats yt 209 if @tpe@ « sede pirat a 
aoe wg | int A ean ratetous To hibes bavgrge) a4 /yoeryace 
-steTaeth Gheinewese eft .oftialenh coon @t biwew aonyqissdeb gommste 
tebenaana hi) beoelanl aou iV ont Geel vofsous al iv 
waetpa to, eteeeol vs rw i< sit Gin’ Yeaped aunnyeg az ae 
io lea ? 2 ee fi OOb PVG lap eves. tlw el eben por tatal al 
mans 5) O0¢RtE jit Otphar var, @otpaat esl) wilt 9 Pasge ; 
od 14 esndlvws ob (pele Ghia qpiumédt Dib sad? oA: A ae a 


eeewe. inde .ali «a eallbiegl soeegaedga th hive exmnriquot aly 


oli oo) tirey. .@l al leve yor « =i sulan dandl si bt sariok 06 - — 


734 2u6y mats fal: Wail al Ao7 OG Aint? tt dorg Zep suhgHNh T= 1 Ane mite ; 
YOtm) €h: (ore ‘tay. ied Vi yey eouit ai onan | . ide io 
Asam ep ne5. ont Hi VATb uaa tay oldman Prd on CAT 
hy 


blues. 2301Ke amt ..2oniertils, To Gyenta ore appa 
4 

en Liog a uate, aap eigen 518i 

ebegn AUTheen er 4 %e 


AG aap ete 


90 


lie in the extension of the theory@tonthe various ‘real-life’ situations 
discussed in Chapter VI. This would also assist the experimentalist 

in testing the different alternatives, The behaviour of the turbulence 
Statistics is an important question to be answered. The hypothesis 
relating shear stress to turbulent energy is in need of verification, 


for it is a powerful tool in studying flow development. 
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APPENDIX B 
DETAILS OF NUMERICAL SCHEMES OF P, A, TAYLOR (1972) 


The equations of the mixing -length, Glushko, and modified 
Peterson models are all essentially parabolic in nature and can be solved 
with a finite difference representation moving forward in x one step at 
a time, The mixing-length model gives rise to a single parabolic second 
order partial differential equation, essentially in U, to be solved 
simultaneously with the continuity equation, which plays a secondary 
role in this context, The Glushko model gives two simultaneous parabolic 
partial differential equations together with continuity, The modified 
Peterson model involves a single parabolic partial differential equation, 
the energy equation, to be solved simultaneously with the two first order 
equations, those of momentum and continuity. To avoid discretization 
errors, the transformation to a vertical coordinate ”=1n (242z,,)/2, 


is employed, All speeds are scaled with respect to u,, and lengths 


1 


with respect to z For computational stability the 'diffusion' terms , 


9° 
that is, oe are represented implicitly, A system of equations 

can then be set up with only a tridiagonal matrix to invert. The finite 
difference grid is shown in Fig, B-l. The grid spacings H and STP as well 
as AG are changed in a block structure as the integration proceeds 


downstream, as shown in Fig, B-2, 


B,1 The Mixing-length Model 


The horizontal momentum equation combined with the mixing length 


hypothesis in JS coordinates is 
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Figure B-1, The finite difference grid. 


Figure B-2. Computation blocks, showing change in H and gradual 


increase ini, 
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1 j i= 
This equation represents N simultaneous linear equations in UN., the 
value of U at the downstream step, assuming that U and W are known, 


In matrix form 


Do 

[ 1 2 eS By 
D. OD 
OD. D O 9 UN, By 
OD. Deo Ds UN re B. 
, j j 
D 

oD, Dy UN By 


The diagonal and off-diagonal elements of the tridiagonal matrix are 
given by 
D,=U,/STP 4+ 2K-M DU. fH 
J A df 
Ob = ew 
J peed 
The elements of matrix B are given by 
B.-U*/sTP - M.W.DU, + k°M.DU.D2U, 
J J Joe! Bilin tal J 
The lower boundary condition makes no changes, but for j=N the outer 


boundary condition requires that 


2 Z 
D =U, /STP +k My DU, /H 


N 
4 Bs 2 
= = D D2 ALF /H 
B= Uy/STP MW DU, +k M U6 Unt LF /H’) 
ae = (US UN = UN. . OD _ is unchaneed, 
where ALF =Uy 4 N Nel N N ang 


The tridiagonal set of linear equations can then be solved using a 
direct 'two-pass' Gauss elimination procedure, With the values for UN, 
at X+STP the values of W are computed using backward differences for 
-1 -1 
—=W, _-(H/2(M, (UN.- U,)+ M, , (UN, ,-U. _))/STP 
‘ie puns oe j > jee jek 
This method was found by Taylor to be stable for all values of STP, 
The values of the control parameters and step and grid sizes used in 


this application are in Table B-1, Computation took about 20 s of CPU 


time on the University's IBM 360-67, 
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B,2 The Glushko Model 
With the coordinate transformation the governing equations can 


be written as 
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In finite. difference form, with a as the value of the energy at X+STP, 
the first and third equations become 
gba’ eas + M.W.DU. 
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For the momentum oo the matrix elements are given by 
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The lower cee condition does not affect these expressions, but the 


outer boundary condition leads oo8 
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with OD, unchanged from the above, 
For the turbulent GUESS equation the matrix elements are given by 
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The boundary condition at the ground, = = 0 , can be imposed in the 


form EQ = Fy EN ee. This necessitates a change in Djs namely, 


D, =U, /StP+ (1/2)k M gil? py? 
1 Pay 
OD, and By are calculated as before, but with Eo Fy The outer boundary 
condition is imposed as Eu p= 1.0/a. This leaves De and OD, unchanged, 


but By! has an additional tern, 


By = By (as before) + (1/2) ME Pag 


Both systems are solved separately for oN and se and u, is computed as 
in the mixing-length model, 

Taylor found that the system was not stable for all step sizes, 
and shorter steps were needed than in the mixing-length case, The values 
used in this application are given in Table B-2, The computations 


required 300 s of CPU time on the IBM 360-67, 


B,3 The Modified Peterson Model 


In J coordinates the governing equations can be written 
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In finite difference form with a as the stress at X+STP the turbulent 


energy equation becomes 
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B.=U,T,/STP-M.W,DT,+ a TM,W,-a TM Me 
J tal J aq m a es J J 
+(1/2)u «crt! D2T + pty?) 
a) J J J J 
The boundary condition at the ground, ®t -~o6 , is imposed as Ty =T)> 
no) 
= UN 
TN5 ? so nay 
= S 
D,=U,/STP + (1/2)kM,T /H 
The outer boundary condition, Tit 1,0, requires that Be be modi fied, 
B_=B (as before) + (1/2)kM ghl2 py? 
N N NN 


With the values for ae the speed uM, can be computed from the momentum 


equation according to 
UN,.=U,+ (1/2)M,STP(DT,+DIN,)/U. -STP M.W.DU./U, 
i J J J J J fie Pa aa | 


‘tos TN.) /2H. 


where DIN, = (TN. 
J a 
Once again Taylor found that the scheme was not unconditionally 
stable, smaller step sizes being needed than either the Glushko or 
mixing-length models, The values used in this application are given in 
Table B-3, The computations required about 450 s of CPU time on the 


IBM 360-67, 
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APPENDIX C 


NOTATION 


constant of proportionality between the standard deviation 
of the vertical wind and the friction velocity 

also the area of a lot containing roughness elements 
constant of proportionality between the shear stress and 

the mean turbulent energy 

arbitrary constants in various expressions 

constant of proportionality between the standard deviation 
of the longitudinal wind component and the friction velocity 


drag coefficient relating the wind speed to the friction 
velocity 


drag coefficient for the first surface defined in terms of 
the wind at 8m 


drag coefficient for the second surface defined in terms 
of the wind at 1/2 m 


specific heat of air at constant pressure 

zero-plane displacement 

also the depth of the modified layer in the Panofsky- 
Townsend theory 

base of natural logarithms 

mean specific turbulent kinetic energy 

fluctuating part of the turbulent energy 

exponential integral, defined at first occurrence 
blending function for horizontal velocity 

blending function for shear stress 

blending function for vertical velocity 

universal function of stability for vertical wind 
fluctuations 

universal function of stability for longitudinal wind 
fluctuations 

universal functions of the parameter % in the Townsend 
similarity theory 

acceleration due to gravity 

functions of the roughness ratio in the internal boundary 
layer growth equation proposed by Shir 


functions used in the iteration to obtain a least-squares 
solution for the profile parameters 


heat flux 

also channel height in Jacobs! work 
height of the internal boundary layer 
height of the new equilibrium layer 
crop height 
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incomplete gamma function used in Dyer theory 
indices 

eddy diffusivity 

exchange coefficient for heat 

exchange coefficient for momentum 


exchange coefficient for turbulent energy 
von Karman's constant 


constant of proportionality between mean scale of turbulence 
and height 


constant of proportionality between dissipation length and 
height in Glushko model 


mixing-length 
length scale in Townsend similarity theory 
Monin-Obukhov length 


related length (= Ky/Ky L) 
mean scale of turbulence in Glushko model 
dissipation length 


roughness change ratio (= Z91/Z92) 

roughness change parameter (=1nm 

number of measurement levels 

also number of roughness elements in a specified area 
index for iterations 

also exponent in one hypothesis of R.J. Taylor theory 
parameter in Townsend similarity theory 

function of height only occurring in Jacobs! combined 
equation of motion and continuity 

Richardson number 

parameter in Businger-Dyer-Paulson diabatic profile 
as a subscript, indicates reference level 

specific or lot area of roughness elements 

silhouette or cross-sectional area of typical roughness 
element 

time 

mean temperature over the profile 

fluctuating part of temperature 

mean horizontal wind 

fluctuating part of horizontal wind 

mean horizontal wind over first surface 

mean horizontal wind over second surface 

friction velocity at surface 

triction velocitva au LiAnstweuniace 

friction velocity at second surface 

velocity scale in Townsend similarity theory 
dimensionless horizontal velocity (= u/ux)) 
fluctuating part of lateral wind 

mean vertical wind 

fluctuating part of vertical wind 
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mean vertical wind over first surface 

mean vertical wind over second surface 
dimensionless vertical wind (= w/u.]) 

downwind distance from roughness discontinuity 
dimensionless downwind distance (= x/z 9) 
vertical coordinate, height 

roughness of first surface 

roughness of second surface 

height of surface boundary layer 

dimensionless height (= z2/z,9) 


symbol introduced to shorten the expressions used in the 
iterative scheme for the profile parameters 


empirical constant in diabatic profile theory 
streamline displacement 

dissipation rate 

transformed vertical coordinate for numerical computations 
(= ine z+2 5) / 249) 

parameter in Townsend similarity theory (= z/1,) 
potential temperature 

also concentration of diffusing entity in Dyer theory 
mixing-length at outer edge of boundary layer 
viscosity 

dummy variable 

density of air 

standard deviation of longitudinal wind component 
standard deviation of lateral wind component 

standard deviation of vertical wind component 
kinematic shear stress 

kinematic shear stress over first surface 

kinematic shear stress over second surface 


value of shear stress at first surface (=u) 9 
value of shear stress at second surface (=uxg9) 
non-dimensional wind shear 


also flux in Dyer theory 

equilibrium flux value in Dyer theory 

centroid of region of streamline displacement 
streamfunction 

value of streamfunction along the streamline CD 
two-dimensio:.1l vorticity 
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